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ABSTRACT 

A  preliminary  design  of  a  coal  burning  25,000 
H?  naval  ship  propulsion  plant  was  completed.   The  system  con- 
sists of  a  combined  Closed  Brayton  cycle  and  bottoming  RarJcine 
cycle.   Coal  combustion  capability  is  achieved  by  using  a 
Pressurised  Fluidized  Bed  (PFB)  heater.   A  thermodynamic 
analysis  and  volume  estimate  was  conducted  for  various  Brayton 
cycle  compressor  pressure  ratios  (  7TC  )  . 

The  results  show  that  a  maximum  thermodynamic 
efficiency  (  ff  )  is  achieved  in  the  range  of  TTC  =  5  to  6 . 
Total  volume  on  the  other  hand  has  a  minimum  around  T7C  =  3.5 
and  climbs  steeply  above  that.   A  trade-off  exists  between 
efficiency  and  volume  with  a  good  design  range  of  7TC  =  U>   to  5  • 
The  design  point  was  ultimately  chosen  at  77"^=  4-. 5  with  a  Of-  3 
AQ.  ' 

The  design  of  a  modern  coal  burning  combined  cycle 
offering  high  efficiency  at  both  rated  and  part  load  is  within 
the  reach  of  present  technology.   Heat  exchangers  and  turbo- 
machinery  design  are  well  within  the  state  of  the  art  and  PFB 
technology  is  expected  to  come  on  line  within  the  next  few  years. 
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ABSTRACT 

A  preliminary  design  of  a  coal  burning  25 , 000 
HP  naval  ship  propulsion  plant  was  completed.   The  system  con- 
sists of  a  combined  Closed  Brayton  cycle  and  bottoming  Rankine 
cycle.   Coal  combustion  capability  is  achieved  by  using  a 
Pressurized  Fluidized  Bed  (PFE)  heater.   A  thermodynamic 
analysis  and  volume  estimate  was  conducted  for  various  Brayton 
cycle  compressor  pressure  ratios  (  7TC  )  . 

The  results  show  that  a  maximum  thermodynamic 
efficiency  (  Vf  )  is  achieved  in  the  range  of  7TC  =  5  to  6. 
Total  volume  on  the  other  hand  has  a  minimum  around  TTC  =  3.5 
and  climbs  steeply  above  that.   A  trade-off  exists  between 
efficiency  and  volume  with  a  good  design  range  of  7TC  =  ^  to  5« 
The  design  point  was  ultimately  chosen  at  77c  =  ^.5  with  a  Of  = 
AQ.  ' 

The  design  of  a  modern  coal  burning  combined  cycle 
offering  high  efficiency  at  both  rated  and  part  load  is  within 
the  reach  of  present  technology.   Heat  exchangers  and  turbo- 
machinery  design  are  well  within  the  state  of  the  art  and  PFB 
technology  is  expected  to  come  on  line  within  the  next  few  years. 
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I.   INTRODUCTION 

1.   Purpose*   Increasing  prices  and  decreasing  supplies  of 
fuel  oil  has  created  a  need  for  highly  efficient  ship  pro- 
pulsion systems  capable  of  utilizing  fuels  other  than  oil. 
The  purpose  of  this  thesis  is  to  perform  a  preliminary  design 
of  a  combined  cycle  propulsion  plant  equipped  with  a  multi- 
fuels  combustor  capable  of  burning  coal  or  oil. 


2.  System  Description.   The  system  investigated  is  a  closed 
regenerative  Brayton  Cycle  combined  with  a  bottoming  Rank in e 
cycle.   Working  fluids  for  the  cycles  are  air  and  water 
respectively.   The  heater  for  the  Brayton  cycle  is  a  pressur- 
ized fluidized  bed  (PFB).   The  heater  cycle  is  essentially  an 
open  regenerative  Brayton  cycle.   This  thesis  was  conducted 
considering  the  heater  only  in  its  coal  burning  mode  of  opera- 
tion. Work  done  by  the  Brayton  and  Rankine  cycles  will  drive 
a  common  shaft  through  a  double  reduction  type  gear.   An  enthalpy- 
enthropy  diagram  is  shown  on  Figure  1 . 


3.  Scooe. 

a.   Objectives.   The  objectives  of  this  thesis  are  to: 
(1)   Perform  a  thermodynamic  analysis  of  the  system 
for  various  closed  Brayton  cycle  pressure  ratios  utilizing 


h  f 


Figure  1 . 


Closed  Brayton  Cycle 


TT> 


realistic  values  for  component  efficiencies  and  losses. 

(2)  Estimate  major  component  volumes  and  thus  obtain 
a  total  volume  for  set  of  parameters  listed  in  (1)  above. 

(3)  Select  a  set  of  parameters  for  design  based 
upon  the  above  results. 

(*0   Discuss  the  technical  feasibility  and  implications 
of  the  design  as  a  means  for  naval  ship  propulsion. 

b.  Design  Criteria. 

(1)  The  system  parameters  selected  for  design  were 
those  which  gave  the  lowest  total  volume  within  an  acceptable 
range  of  thermal  efficiency. 

(2)  An  acceptable  range  of  thermal  efficiency  is 
defined  as  1  1%   of  the  maximum  efficiency. 

c .  Limitations  to  Scope. 

(1)  No  component  or  system  testing  was  carried  out. 

(2)  The  system  was  designed  with  the  aid  of  a 
computer  utilizing  performance  parameters  of  currently  avail- 
able equipment. 

(3)  The  concept  of  pressurized  fluidized  bed  (PFB) 
combustion  has  not  been  proven  in  a  shipboard  environment. 
The  thesis  was  conducted  based  on  the  assumption  that  a  PFB 
system  would  be  put  to  sea  in  the  near  future. 
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II.   THERMODYNAMIC  ANALYSIS 

1 .  The  Brayton  Cycle 

a.  Working  cycle.  This  is  a  closed  Brayton  cycle 
utilizing  air  as  the  working  fluid.  Air  was  chosen  since 
this  is  proposed  for  surface  ship  application  and  would  thus 
eliminate  the  need  for  carrying  a  reserve  supply  of  gas  if 
something  like  Helium  or  C02  were  used. 

Advantages  of  a  closed  cycle  over  an  open  cycle  are 
the  following: 

(1)  Compressor  suction  pressure  is  no  longer  limited 
to  ambient.   By  increasing  the  working  pressure  of  the  cycle 
a  significant  reduction  in  engine  size  can  "be  achieved  for  a 
specified  output. 

(2)  Higher  working  pressures  also  increases  the  heat 
transfer  capacity  of  the  working  fluid  thus  resulting  in 
smaller,  more  effective  heat  exchangers. 

(3)  -he  use  of  a  high  pressure  fluid  gives  the  added 
advantage  of  regulating  the  cycle  pressure  level  by  varying 
the  inventory  of  gas  in  the  loop.   If  this  is  done  in  a  manner 
which  kept  the  maximum  temperature  and  speed  constant,  the 
turbomachinery  vector  diagrams  would  remain  essentially  un- 
changed.  The  result  is  little  or  no  change  in  compressor  or 
turbine  efficiencies  with  varying  load,  thus  maintaining  high 
part  load  efficiency. 

A  Temperature  -  entropy  (T-s)  diagram  of  the  working 
cycle  is  shown  in  Figure  2. 


T  * 
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Figure  2. 

compressor  suction,  cooler  discharge 
compressor  discharge,  regen  inlet  (cold) 
regen  outlet  (cold),  heater  inlet 
heater  outlet,  turbine  inlet 
turbine  outlet,  regen  inlet  (hot) 
regen  outlet  (hot),  WHB  inlet 
VJHB  outlet,  cooler  inlet 


Qs:    heat  input  from  heater  cycle 

■ 

QR:   heat  released  to  Rankine  cycle 

« 

Qc:   heat  released  to  sea  water  cooler 
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The  analysis  was  conducted  by  assuming  values  for 
T- ,  TV,  7cfp)(comPressor  polytropic  efficiency),  tj^   (P)  (turbine 
polytropic  efficiency),  ER  (regenerator  effectiveness),  and 
[p~jT    (total  pressure  drop). 

The  specific  heat  was  assumed  to  vary  with  temperature 
according  to  the  following  relation: 

Cp(T)  =  .2^75  -  (3-759xl0"5)  T  +  (5.lc6xlO"8)T2-( 1 .231X10"11 )T3 

for     T  in  °R 

Cp  in  BTU/lbm-°R 

The  following  are  defined; 

TT^   -   compressor  pressure  ratio  (P^/P.,) 

jft   -   turbine  pressure  ratio  (P«-/Pj«) 
(4rL        -       f°   total  pressure  drop  around  closed  Bray  ton  loop 

* 

WB   -   power  (net)  from  closed  Brayton  cycle 

m    -   mass  flow  of  working  fluid  in  Brayton  cycle 

T.    -   temp  in  °R  at  point  i 

h.    -   enthalpy  at  point  i 

R    -   gas  constant 

P.    -   pressure  abs .  at  point  i 

V    -   specific  volume 

ER   -   regenerator  effectiveness 

The  basic  procedure  was  to  determine  T-,  T~,  T.,  Tr, 
^inp  and  WB  "for  var^-ous  pressure  ratios  (  TT*  )  . 
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Determination  of  T2  (ref.  Dixon): 


Small  stage  polytropic  efficiency  is 

W>  dh  CPdT 


since  for  an  isentropic  process  T,   =  0=dh.  -  vdp 
Substituting     IT  -   RT/p  gives 


Cp  '  p  *  dT 


thus, 


T, 


It   =  :r^- 


5 


-     ^^   ?    ?c(p)   I?, 


By  substituting  the  equation  for  C   =  C _(T)  one  gets 
.2475  ln/T2\    3.759xlO"5(T2-T1)  +  (2 .553xlO'8) (T22-T12) 

(4.103X10"12)  (T23  -  Tt3) 


=  R     /. 


H?) 


M   TTc 


T2  can  then  be  determined  iteratively 


Determination  of  T,i 

In  a  manner  similar  to  the  compressor  above, 


*  c 


^T  =^^(f) 
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After  integrating,  this  "becomes 

.2475  ln/T5\-   (3.759xlO'5)(T5-T^)+(2.553xlO"8)(T52-T42) 
(4.103xl0"12)    (T53   -  T^3) 

Tc   can  then  be  determined  iteratively. 

Determination  of  T-  and  T/-: 

T~  and  T^  can  be  found  from  an  assumed  regenerator 
effectiveness  (ER)  using  the  following  relation: 


ER  =  T5  -  T6  .  T3  '  T2 
l5         2     x5    2 
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b.   Heater  Cycle*   This  is  an  open  Brayton  cycle  also  using 
combustion  air  as  the  working  fluid.   The  heater  itself  is 
a  pressurized  fluidized  bed  utilizing  coal  as  fuel.   A  T-s 
diagram  of  the  heater  cycle  is  shown  in  Figure  3» 

Major  components  of  the  cycle  consist  of  a  heater,  re- 
cuperator and  turbine  driven  compressor.   The  turbine  is 
powered  by  exhaust  gases  leaving  the  recuperator. 

The  analysis  was  conducted  by  assuming  values  of  T,.t 

T7h'  ?£-(p)»  !?-f(p)'  and  0  (equivalence  ratio).   Both  the 
exhaust  products  and  combustion  air  were  assumed  to  be  perfect 
gases  with  specific  heats  varying  in  accordance  with  the  re- 
lation given  in  l.(a). 

The  following  are  defined; 

Vcu        -       heater  compressor  pressure  ratio  Pph^lh^ 
TTfj,    -   heater  turbine  pressure  ratio  (^v/P^vJ 


m, 


h 


total  pressure  drop  around  heater  loop 


m  -  mass  flow  of  air  in  heater  cycle 

a  " 

m-  -  mass  flow  of  fuel  in  heater  cycle 

(F/A)s  -  stoichiometric  fuel  -  air  ratio  of  fuel  used 

LHV  -  Lower  Heating  Value  of  fuel  used 

Tih  "  absolute  temp  at  point  i, 

h.,  -  enthalpy  at  point  i, 

/  H  -  efficiency  of  the  heater  cycle 

EREr  -  recuperator  effectiveness 


Ik 


T  /T 


ambient ..   lh 


QIN  =  mf  LHV 


Figure  3« 

point  lh:  compressor  inlet 

2h:  compressor  outlet,  recup  inlet  (cold) 

3h:  recup  outlet  (cold),  heater  inlet 

^h:  max  air  temp  in  heater 

5h:  heater  outlet,  recup  inlet  (hot) 

6h:  recup  outlet  (hot),  turbine  inlet 

7ht  turbine  exhaust 


Q^:   heat  released  to  working  cycle 


Qin:   heat  in  =  mf  LHV 
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Assume  the  heater  cycle  to  be  a  control  volume  (C.V.) 
as  shown  in  Figure  ^. 


mfLHV 


> 


C.V. 


r 


[ 


*Q 


B 


a  In 


(ma  +  mf)  hyh 


Figure  k. 

QB  leaving  the  C.V.  must  be  equal  to  QB  entering  the 
working  cycle  thus, 


QB  =  mB  (h^  -  h3)  =(ma  +  mf)(hi+h  -  h^) 


Performing  a  heat  balance  on  the  C.V.  gives  us 


QB  ♦  (ma  ♦  mf)  h?h  =  if(LHV)  ♦  ^ahlh_ 


Therefore, 


~[(i +('/>);  0)h7>    -U]/[(%)jLHV] 
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From  a  power  balance  between  the  turbine  and  compressor 


•     t 


WQ  =  WT  (no  net  work) 


raa     5pc    (Th2    -Thl>    "    ^a   +   if)    CoT    <T6h   "  T7h> 


but 


T*j,    *7m,  (ffch) 


zPc/?Jf>) 


T^-uK)^^ 


thus, 


%a(t'-^9ri.'77r«0 


or, 


C,c 77K  [(^)W?"- 1]  ♦  [  /  +  (%),  0] CFf  T7K  = 


V<%)^]CnT%[,-(^)r,]^JT^^ 


1? 


therefore, 


i  - 


p'rh 


Cf}+ 


For  minimum  volume  heat  exchangers  we  want  the  largest 

(Ap/\ 

*~/?/t'      compatible  with   the   above   eqn. 


differentiating , 


I      5  77^ 


K    .1 


aP\ 


Setting 


\    P/Th 


TT 


=    0  and   solving  for  TTch  gives 


ch 


TT 


C-K 


[i  <FA),  *  ]  jfc  ft  flk  r7A  -  c^  t-k)  7  ^/v 


IS 


By  substitution 

Thus   ^Tck  and  (~~f~)     can  be  found  through  known  or  assumed 
quantities. 

Determination  of  T^h>  T^ 

_         CKfcrl  -Lk)  Cc(Tjk  ~T2k) 

where  Ck  =  Crk  (  ^c<  +  ^-f ) 

Cc    =   ^c   "V 
therefore , 

T2h  and  IV,    can  be   found   once  TTC^  is   determined.      Assume 


that 


Tc.     =   T-   +    100°R 
5h  3 
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2.  The  Rankine  Cycle.   This  is  a  bottoming  steam  cycle  that 
recovers  heat  from  the  closed  Braytcn  cycle  that  would  other- 
wise be  discharged  overboard  through  the  cooler.   The  rela- 
tively hot  gas  leaving  the  regenerator  is  pre-cooled  in  the 
Waste  Heat  Boiler  (WHB)  before  being  sent  to  the  cooler  which 
brings  the  temperature  down  to  compressor  inlet  conditions. 
A  T-s  diagram  of  the  Rankine  cycle  is  shown  in  Figure  5« 


Th 


Figure  5« 
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Major  components  of  the  cycle  are  the  WHB,  condenser, 
steam  turbine,  and  feed  pump.  It  is  assumed  that  the  feed 
pump  is  driven  with  power  from  the  turbine. 

The  analysis  was  conducted  by  assuming  values  of  T^g, 
A  T ,  OpO5)  ( pump  isentropic  efficiency),  D-fO5)  (turbine  isen- 
tropic  efficiency),  and  7  WHB  (waste  heat  boiler  efficiency). 

The  following  are  defined; 

h.    -   enthalpy  at  point  is 

T.     -   absolute  temp  at  point  is 

PP    -   pinch  point  on  Brayton  side 
SF    -   saturation  point  on  steam  side 
/5    _   mass  flow  ratio  (m^/nO 

x    -   quality  of  steam  exiting  turbine 
P.     -   absolute  pressure  at  point  is 

tjvr    -   power  (net)  from  the  steam  cycle 

/f    -   thermal  efficiency  of  the  combined  cycle 
nu    -   mass  flow  of  steam  in  Rankine  cycle 

Rankine  cycle  parameters  were  chosen  for  maximum  worlc. 
The  procedure  was  as  follows: 

(1)  Compute  T-g  from  the  relation 

T33  =  T6  -AT 

(2)  Assume  a  value  for  T   .    Tx  <  T   <  662. 0°F 
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(3)  Compute  Top  from  the  relation 

TSF  =  Tpp  "AT 

(4)  Obtain  PSAT »  hs„,  h-c  from  steam  tables. 

(5)  Compute  hji,s  from  the  relation 

h^s=h-s-  '^-f  (  is)  (h^s-h^s  (is)  )   get  J"^s(,s  )  from  steam  tables 

(6)  Obtain  x  (turbine  exit),  h--.  from  steam  tables 

X  o 

(7)  Compute  h2s  from  the  relation 

L   -  L   *  k"Ci$)   "*« 

(8)  Compute  /£  from  the  relation 


!k& 


(9)   Compute  W„  from  the  relation 

*(10)   Iterate  with  T    (step  2)  until  WR  is  maximum. 

(11)  Compute  T7   from  the  relation 

CpfT7).T7  =  Cf(rPP)-TeP  -fiftsp  -\iu)fo 

(12)  Calculate  jjjj   from  the  relation 


y+  =  wB  ♦  wR 

«B 
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III.   VOLUME  ANALYSIS 

1.  Turbomachinery .   An  analytical  procedure  for  the  minimum 
volume  design  of  turbomachinery  was  developed  by  Lee  (ref . 
10).  The  following  equations  and  mathematical  expressions 
summarize  that  portion  of  the  procedure  which  are  applicable 
to  compressor  and  turbine  designs  in  this  thesis.   Volume 
obtained  by  Lee's  method  neglects  the  volume  of  casing,  inlet 
and  outlet  plenums,  and  bearing  assemblys. 

a.   Turbine.   Free  vortex  blading  and  axial  exit  stages 
were  assumed  throughout  the  analysis. 

The  following  are  defined: 

S.F.    -  factor  of  safety 

U+  -  max.  allowable  tip  speed 

0.  -  flow  coefficient  at  tip 

y~  -  hub-tip  ratio  at  inlet 

V/3  -  hub-tip  ratio  at  exit 

C-j.  -  clearance  factor 

V*  -  specific  volume  at  inlet 

Tt  -  max  tip  radius 

h  -  polytropic  coefficient 

K  -  isentropic  exponent 

P»\  -  density  of  blade  material 

°<  -  stress  concentration  factor  at  root 

a  -  amplification  factor 

9  -  correction  factor  for  taoer 
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xi    ~   dist.  from  tip  to  1st  principle  axis 

chord 

"^n  "   normalized  moment  of  inertia 
Vol..  -   volume  of  turbine 


The  volume  was  computed  by  assuming  values  for  SF,  U^,  0+» 
y*  ,  CL ,  K ,  /?m  ,  *  ,  a ,  x1 ,  Ksl  ,  and  0 . 

r, ,  n,  and  y£    are  computed  from  the  following  relations: 


n  - 


*  -fa-/)fc<W 


n 


Compute  C- ,  Col  C~,  C^  from  the  following: 

Cx-  z(A-£T^)/(SF-«) 

C3  =  Z3TK(l-^)""/(iF^) 

c4  =  -ep^a^/z 

where  A  and  B  are  material  dependent  constants.   The  following 
can  then  be  calculated: 
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The  axial  chord  for  a  stage  is  given  by  bL •   For  a 
large  number  of  stages,  summation  of  b,  for  the  whole  turbine 
can  be  approximated  by  an  integral: 


where 


This  integral  can  be  evaluated  numerically  thus, 


Volf  =  27rr/ct  IT 

b.   Compressor.   An  approach  similar  to  that  of  the 
turbine  was  used  for  the  axial  compressor.   Axial  exit  stages 
were  assumed  in  the  analysis. 

The  volume  was  computed  by  assuming  values  for  SF,  U+, 
^iHc,  CL,  K,  D,  p<  ,  a,  xl ,  KI(  ,  9,  and  T  (blade  material 
yield  stress)  . 
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r+f  n,  and  y£    are  computed  from  the  following  relations: 


n 


n  = 


K 


K  -(K-O/Vcfo 


2.  I  /jfc 


,1 


Compute  C1  ,  C-*  and  C~  from  the  following: 


C3   =  2  Vr  T  *V 

Then  calculate  the  following: 

N  -  ct[(i  -y)/0+»)](fS  itf^/ifo*  t.rzisy1-) 


D 


--G^O-^j 


j>j)  =  -z  c2  y 
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As  with  the   turbine   define 


'J 


where 

This  integral  can  be  evaluated  numerically  thus 

voi  =  zTrn,zcL  IT 


c 


2.  Heat  Exchangers.   The  volumes  of  the  Regenerator,  Re- 
cuperator, Cooler,  and  Heater  were  computed  using  the  method 
of  minimum  volume  design  developed  by  Lee  (ref.  10).   In 
essence  Lee  found  that  the  optimal  diameter  ratio  (D./D)  for 
these  heat  exchangers  can  be  established  as  a  simple  relation 
to  the  pressure  ratio  {Tf    )  of  the  cycle.   The  results  of  the 
optimally  distributed  pressure  drops  (  -=-  )  as  reported  by 
Lee  were  also  studied.   Parametric  equations  relating  pressure 
drop  to  pressure  ratio  were  developed  and  used  in  this  analysis 
No  consideration  has  been  made  on  computing  the  volume  of  the 
shell,  headers,  or  other  appendages.   Lee's  relations  were 
developed  from  basic  heat  transfer  and  pressure  drop  relation- 
ships. 
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a.   Regenerator*  The  regenerator  is  assumed  to  be  a 

counter  flow  shell  and  tube  heat  exchanger.   High  pressure 

is  inside  the  tubes  as  is  the  case  for  current  design  prac< 
tices. 


The  following  are  defined: 

f4r)    -   $  pressure  drop  across  the  regenerator 
•P/R,!i      (both  steams) 


m 


Reg 


diameter  ratio  for  regenerator 


/u  -  enthalphy  difference  across  the  interior 

'Ks-g  stream  of  Reg. 

D  -  tube  inside  diameter 

AT  -  log  mean  temp,  difference  of  hot  and  cold 
streams 

P„  -  Prandlt  number 
r 

jp  -  density 

/<  -  dynamic  viscosity 

k  -  thermal  conductivity 


/AP\        /DK\ 
Values  for  (TyRe-  and  [TJ~j  oexT  can  be  estimated  from 

the  following  relations: 


(^Reg  =  -135  -  •°28^c  (for(fJr  =  •*•  ER  =  -88) 


(-5s)  Reg  "  -557  ♦  -0933  7T, 
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Subscript  1  refers  to  quantities  inside  the  tubes 
whereas  subscript  2  refers  to  quantities  outside  the  tubes. 


A„    can  be  computed  as  follows: 
Reg 


A*.,  a(%)  ^  Yk  ^e  ^«  '"+l  [  '  +(*)» 


-,  2.4/ 


where 


r*l7.Z*kt      A    ft.    /? 

y  -  n      s>     at 


The  volume  can  then  be  computed  as  follows: 

/AP\    -0.4.1 
V01Reg  "  AReg  fT/Reg 


b.   Recuperator,   The  recuperator  is  similar  in  design 
to  the  regenerator  only  it  refers  to  the  heater  cycle. 


The  following  are  defined; 
!v  p  /Rec 


d  y Rec 

/AP\ 


AP\ 

7/PFB 


AK 


1  Rec 


$  pressure  drop  across  the  recuperator 
(both  streams) 

diameter  ratio  for  recuperator 

%   pressure  drop  across  the  heater  side 
(outside  tubes) 

%   pressure  drop  across  the  PFB 

Enthalpy  difference  across  the  interior 
stream  of  recup. 


29 


Values  for(""p"/Rec  and("5L)Rec  can  be  estimated  from  the 
following  relations 


=  /c 


Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
the  tubes  respectively.   AR   is  computed  as  follows: 

where  V  and  J  are  as  previously  defined. 
Then  compute  the  volume 

M?\    -0.41 


VolRec   kRec{) 


c.   Cooler.   The  cooler  is  assumed  to  be  a  shell  and  tube 
heat  exchanger  with  the  gas  inside  the  tubes  and  salt  water 
outside. 


The  following  are  defined: 

/—) 

l^P/Clr   -   %   pressure  drop  across  the  cooler  (air  side) 

]p~/Clr   -   diameter  ratio  for  cooler 

^K,  Clr   -   enthalpy  difference  across  the  interior 
stream  of  cooler 

Y 
WKB2    -   sp.  vol.  of  working  fluid  at  WHB  exit 
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/A?\ 

l  p  /WHB     -   #  pressure  drop  across  the  WHB  (air  side) 

fp~)i       -   ^  pressure  drop  across  the  heater 
^     '         (inside  tubes) 

^pOl)  -       polytropic  cooling  water  pump  efficiency 

V;.      -   sp.  vol.  of  cooling  water 

Values  f or  \r  /Clr  and\vj}/Clr  are  estimated  from  the 
following  relations: 

/JP)        fc/i£\   .  /AP\         JAP)        J&?\ 


Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
the  tubes  respectively.  Ap,   is  computed  as  follows: 

where  v  ,  Y   are  as  previously  defined  and 

c  =  (7«,)(^^"l"(%,)(%) 
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The   volume   is   thus 

VOW  "Aclr  (f  W^ 


d.  Heater,  The  heater  is  assumed  to  be  a  shell  and 
tube  heat  exchanger  with  the  working  gas  in  the  tubes  and 
combustion  gas  out. 

The  following  are  defined: 


J>K 


57  htr 


diameter  ratio  for  heater 


dh  htr   ~   enthaphy  difference  across  the  interior 
1  stream  of  heater 

Values  of(A£\  lAL)..  and  (Dh)     are 

\rMrx,   rp/htr2'        d"  htr 

estimated  from  the  following  relations: 

y^r,     "    I  P  /Wr*   '  (j>  Afr 


/A? 


where 


-1.7S  ,„         N-  zS 


/wv«.r  hs  m(%) 


-  .O^^TTr    +.2/+ 


5  /ht 


Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
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the  tubes  respectively.   Compute  A,,   as  follows: 


Akfr  »wh  ^^^^m^ii^mj'^c 


where  V   ,  X     ,     °^       are  as  previously  defined. 


Then  compute  Volu+.p 


Vhtr  "  Ahtr  (  p  ) 


htr0 


e.   rtaste  Heat  Boiler.   The  boiler  is  assumed  to  be  a 

once  through  cross-counterflow  type,  as  shown  in  Figure  6. 

High  pressure  is  inside  the  tubes  thus  saturation  pressure 

must  be  greater  than  the  working  pressure  of  the  closed 

Brayton  cycle 

gas  ou 

1 


H  H  1  *  t  »  >  |  I 


Ci  i  i  i  ■ 

+  \  I  \   i  I  t  I  *  K 

Nhii'tlii 

CM  1  t  \  1  l  I  li  \r 
1  I  l  I  i  l  I  I  I  IK 

a\  >  >  1 1 1  i  iur 
^n  » i i  i  h  h i 


water  in 


steam  out 


Figure  6. 
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Boiler  volume  is  estimated  using  heat  transfer  relations 
and  volumes  from  similar  U.S.  Navy  COGAS  (Combined  Gas  Turbine 
and  Steam)  designs  (ref.  11). 

The  following  are  defined t 

C  „   -   ave.  specific  heat  on  steam  side 

ps  F 

C     -   ave.  specific  heat  on  air  side 

NTU   -   Number  Transfer  Units 
#Ps   -   number  of  passes 
AHT   -   Heat  transfer  area 

U     -   overall  heat  transfer  coefficient 
V«    -   total  volume 

From  a  DTNSRDC  (David  Taylor  Naval  Ship  Research  and 
Development  Center)  study  (ref. 11)  the  following  relations 
are  obtained: 


,  „    _    CM/A/              ^1«   C?s 

Of- 

mc  Cp 

u            CHflX             m;    CP 

m&  tf>s 

£R*frntec*-,)/(yim- 

0 

NTU  -CHIN 


A*T     = 


U 
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For  similar  designs  assume  A„T  *<   V, 


UI  *  UII 


thus    VT     CMINT     NTUT 


VI];    CMINjj    NTUjj 


vi  °mvii  \  cmint  ntut 

CMINjj  NTUXI 


The  term/VTT         \  can  be  evaluated  from  known  data 

icMTNTT  ru. 


on  previous  designs  to  give 


VWHB  =  ll,53  *  CMIN  *  NTU 


But  it  should  be  noted  that  the  above  equation  assumes 

/A?  \ 
constant  {~J~ J  WHB  and  boiler  power.   If  as  before  Vol  •<■  A„,„ 

then  from  heat  transfer 

AHT  °<  f   «  fi*s  &k 

thus 

Correcting  the  above  equation  for  power  gives 

VWHB  =  ,0017  CMIN  *  NTU  *  ^   *  ^ 
As  stated  previously  Vol  °<  (""£""/ 
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Correcting  again  for  pressure  drop  yields 

VWHB  =  .000^5  CMIN  *  NTU  *■  ms  *-  &hs  (^r~J         (ft3) 

CMIN  in  BTU/s  -  °R 
m  in  lbra/s 

4hg  in  BTU/lbm 


f.  Condenser.   The  condenser  is  assumed  to  be  a  shell 

and  tube  crossfiow  heat  exchanger  with  salt  water  inside  the 
tubes  and  steam  outside. 

The  following  are  defined. 

D  -   inside  tube  diameter 

N  -   number  of  tubes 

L  -   length  of  tubes 

g  -   Gravitational  acceleration  constant 

APo  -   pressure  drop  on  the  salt  water  side 

f  -   friction  factor 

f c  -   density  of  cooling  water 

V  -   velocity  of  cooling  water 


A   -   cross  sect,  free  flow  area  on  water  side 
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w  -  flow  rate  of  cooling  water 
c 

q  -  heat  flow  in  condenser 

C  -  w  (x) specific  heat  of  cooling  water  (C  ) 
c         c  p 

AT  -  temp,  difference  across  the  cooling  water  side 

X  -  steam  quality  entering  condenser 

h~  -  enthapy  difference  between  sat.  liquid  and 


fg 


sat.  steam 


m    -   mass  flow  of  steam  entering  condenser 

9 


From  basic  heat  transfer  relationships 


where 


AHT  =  ^DNL 


H-^ 


Yi\* 


2 


therefore, 


ht     fW- 


where  ur 

A      =  — 
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thus 


^wr 


=   23  ^Pc  ^U   ms  % 


For  similar  designs  assume  V^  •<  A 


HT 


thus 


fas  x)x 


fc  I  for  similar  designs  assuming  same  A   P 


therefore 


V. 


rv\<  *)L 


The  quantity  in  brackets  can  be  evaluated  from  known 
data  on  previous  designs  (ref.  11)  to  give 


Vcond  ■  27-3  V     "■ 


m  in  lbm/s 
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IV.   DESIGN  OF  A  25,000  HP  PLANT 

A  proposed  design  of  a  25,000  HP  coal -fired  combined 

cycle  was  completed.   Coal  combustion  is  accomplished  using 

a  pressurized  fluidized  bed  (PFB).  Design  parameters  were 

eventually  chosen  after  investigating  the  effect  of  changing 

pressure  ratios,  varying  ,— pr J     in  the  Brayton  cycle,  varying 

ER  (regenerator  effectiveness),  different  compressor  inlet  temp. 

(T.,)  and  steam  condensing  temp.  (T-s),  and  variations  in  Brayton 

working  pressure  (WP)  and  saturation  pressure  (P_-   )• 

mm . 

a.   The  Fluidized  Bed.   The  pressurized  fluidized  bed 
operates  at  a  pressure  of  3  "to  10  atmospheres.   Pressure  is 
maintained  by  an  exhaust  driven  compressor.   For  proper  com- 
bustion incoming  air  must  be  preheated  by  the  recuperator  to 
at  least  700°F.   The  bed  material  is  primarily  limestone  which 
serves  to  capture  sulfer  in  the  coal  and  thus  minimize  SO- 
emissions.   A  bed  temperature  of  1700°F  is  generally  used  to 
keep  NO   emissions  at  an  acceptable  level. 

The  fluid  bed  itself  is  about  U,   ft.  deep.   Preheated 
combustion  air  enters  the  bed  via  a  perforated  bed  plate  under- 
neath the  bed.   Pressure  drop  across  the  bed  is  assumed  to  be 
1095. 

The  PFB  is  used  to  heat  the  working  fluid  in  the  closed 
Brayton  cycle.   The  working  fluid  leaves  the  regenerator  and 
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enters  the  tubes  in  the  PFB.   The  working  fluid  is  first 
heated  by  convection  from  the  gas  leaving  the  bed  and  is  then 
ducted  to  tubes  submerged  in  the  fluid  bed  for  final  heating 
to  the  turbine  inlet  temp  (Tj,).  (Figure  7) 


working  fluid  in(T~^ 


combustion 


air  in  (T^h) 


/1s 


0 


/N 


+  ' 


coal  in s 


exhaust 
A 


Fluid  bed 
(1700°F) 


C 


^ 


-*- 


j 


r 


grid  plate 


-        working  fluid  out 


Figure  7 
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b.   Assumptions,   The  computer  program  was  run  with  the 
following  input  values: 


%  (P) 

= 

.88 

r/t(p) 

= 

.91 

ER 

= 

.88 

;K'0 

= 

.80 

J>t<!<J 

= 

.88 

(?), 

= 

.10 

/  MB 

= 

.9^ 

Tl 

= 

5^0.0   °R 

T/, 

= 

2060.0    °R 

T1S 

= 

5^0.0    °R 

At 

= 

+35-0  °R 

wp 

= 

135.0  psi 

p 

3AT 

s 

200.0  psi   (want  at  1 
than  WP) 

Power 

= 

25000  HP 

Tm 

= 

535.0   °R 

T7h 

— 

760   °R    (limited   by  a 
turbine) 

UrX 

= 

.85 

yY^K 

= 

.89 

^ 

= 

.90   (11$  excess  air) 

In  addition  it  was  assumed  that  scoop  injection  could  be 
used  for  both  the  cooler  and  condenser  at  speeds  greater  than 
10  knots,  therefore  power  requirements  for  salt  water  circulating 


pumps  were  excluded  in  the  calculations. 

c.   Results .   The  computer  output  is  provided  in  Appendix 
3.   Note  that  the  saturation  pressure  for  pressure  ratios  of 
2  and  3  were  below  the  200  psia  specification.   This  is  due 
to  the  fact  that  the  low  pressure  stream  leaving  the  regener- 
ator is  at  a  temperature  below  the  saturation  temp,  at  200psia. 
A  plot  of  Iftf  vs.  TT    and  Vol.  ,  ,  vs .  TT    is  shown  in  Figure 
8. 

Note  that  the  volume  increases  sharply  above    c   ~   ^ 
due  largely  to  the  rapid  growth  in  the  waste  heat  boiler. 
Maximum  *7f   occurs  at  around  TT   =  6.   To  investigate  the 
effect  of  (— L  on  0^   and  Vol.  ,  the  program  was  run  for 


f*P 


-)  T  =  .04,  .06,  .08,  .10,  and  .12  between  T[ c   of  4  and  9- 


The  results  plotted  on  Figure  9  show  that  significant 
gains  in  Vr   can  be  realized  with  very  little  increase  in  V- 

for  the  range  of  4-6  pressure  ratio.   The  effect  of  changing 

Zap  ) 
Eu  (regenerator  effectiveness)  in  this  range  with  a  \~f~L   = 

.05  was  further  investigated. 

The  results  listed  in  Table  1  show  no  increase  in  effi- 
ciency as  ER  is  reduced.   Furthermore,  volume  increases  sharply 
as  a  greater  percentage  of  the  work  is  being  done  by  the  Rankine 
cycle  as  ER  decreases. 

Thus  for  greatest  efficiency  with  minimum  volume,  cycle  operation 
should  be  at  the  highest  ER  practicable. 
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Figure   9 
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Table  1 


ER  t  Yolt  (ft3) 


.88  A52  335 

.80  .423  6k6 

.50  .402  7038 


.88  .464  408 

.80  .454  811 

.50  A32  4350 


.88  .474  724 

.80  .464  IO83 

.50  .447  3741 


.88  .475  1050 

.80  .^71  1396 

.50  A56  3508 


.88  .^73  1^50 

.80  .471  1750 

.50  .460  3435 
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To  investigate  the  effect  of  changes  in  T-  (compressor 
inlet  temp.)  runs  were  made  on  k   different  temperatures  and 
the  results  shown  in  Figure  10.   As  expected  efficiency  de- 
creased with  increasing  T.  with  very  little  change  in  the  range 
of  interest  4-4.5  Vin'      Volume  on  the  other  hand  decreased 
slightly  for  TV     less  than  3»5  but  increased  significantly 
above  that.   Thus  for  best  efficiency  and  least  volume  T- 
should  be  kept  as  close  to  cooling  water  temperature  as  possible. 

Finally,  changes  in  working  pressure  of  the  Bray ton  cycle 
were  investigated.  Decreases  in  working  pressure  may  require 
additional  volume  of  Brayton  cycle  components  yet  at  the  same 
time  allow  a  decrease  in  Rankine  cycle  steam  pressure  (recall 
50°R  pressure  differential  between  boiler  streams).  This  will 
in  turn  increase  cycle  efficiency  yet  further  increase  volume. 
These  results  are  born  out  in  Figure  11. 

d.   System  Selection.   In  accordance  with  the  design 
criteria  the  highest  thermal  efficiency  considered  is  0.48? 
(Figure  11).   Thus  it  is  desired  to  select  the  plant  with  min- 
imum volume  at  0-f- =  .482. 

Figure  11  shows  that  only  those  alternatives  with  Brayton 
cycle  working  pressures  of  70psi  and  lOOpsi  meet  the  efficiency 
requirement.   Of  these  two  the  lOOpsi  case  has  minimum  volume 
at  desired  efficiency. 

However,  the  problem  does  not  end  here.   If  one  looks  at 
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Figure   10. 
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the  steam  quality  leaving  the  turbine  at  the  TT   =  5,  WP  = 
100  psi  point  it  is  of  the  order  of  0.84.   This  is  a  little 
lower  than  normally  desired  for  design.   In  fact,  the  steam 
quality  of  all  alternatives  thus  far  considered  within  the 
range  of  interest  is  low.   If  a  steam  quality  of  order  .87  is 
desired  the  current  design  pt.   can  change  in  one  of  3  ways 
as  follows: 

1)  Raise  steam  condensing  temp  (T-s) 

2)  Lower  minimum  saturation  pressure  (Pn*n) 

3)  or  a  combination  of  the  above  two 

Through  trial  and  error  it  was  found  that  case  1  gave  the 
greatest  benefit  with  the  least  penalty.   Computer  runs  were 
made  in  the  4-5  V     range  with  the  steam  condensing  temp 
raised  to  560°R  for  both  the  70  and  lOOpsi  alternatives.   The 
results  show  a  substantial  increase  in  quality,  slight  decrease 
in  efficiency  and  very  little  change  in  total  volume  (Table  2). 


(psi) 
Working  Pressure 

^~c 

n 

X 

Volt  (tth 

70 

*.5 

.480 

.864 

948.3 

70 

5.0 

.482 

.871 

1140.6 

100 

5.0 
Table 

.^79 
2. 

.862 

956.5 

The  70  psi  designs  have  better  steam  quality  and  0  .  than 
the  best  100  psi  case.  In  addition,  at  4.5  TT  the  volume  is 
lower.  Between  the  two  70  psi  cases  above  a  trade-off  exists 
between  Q  .  and  Vol+.   Assume  that  volume  comparable  to  the  case 
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of  TT     -   5,  WP  =  100  psi,  Tls  =  5^0°R  discussed  above  (Figure 
11)  is  desired.   If  this  is  the  case  then  the  new  design  point 
is  at  a  Brayton  cycle  working  pressure  of  70  psi,   IT  =  ^.5t 
T1S  =  56o°R .   The  results  of  this  design  are  summarized  in  the 
following  computer  printout: 
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BRAYTON  CYCLE:   PRESSURE  RATIO  =4.50 
WORKING  PRESSURE  =70.0 
TOTAL  PRESSURE  DROP  =  .05 
MASS  FLOW  (LBM/S)=206.93 


STATE 
1 
2 

3 

4 

5 
6 

PP 

7 

COMPONENT 
REG  EN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(DEG  R) 
540.0 
876.0 

1403-9 
2060.0 

1475.9 
948.0 
836.2 
809.2 
DPP  DHD 
.010  .950 


SUMMARY -.THERMAL  EFF=0.480 

HEATER  EFF=0.933   (PHI=0.90) 

SFC ( LBM/HP -HR ) =0 . 440 

VOL  TOTAL(CU  FT ) =  948.3 

POWER (HP) =25000.0 

WORK   FRAC    BY   STM   CYCLE =0.1 01 


VOL(CU  FT) 

26.5      EFF=0.88 


.009       .490 

.0006 
.030 


5*5 

193.1 

29.7 


EFF(P=0.91 

EFF(P)=0.88 

EFF=0.94 


HEATER    CYCLE:       PRESSURE  -RATI0=4  .  79 
TOTAL   PRESSURE  DROP   =.43 
AIR   FLOW(LBM/S)=    36.85 
FUEL  FLOW(LBM/S)=    3-06 


STATE 

TEMP(DEG   R) 

1H 

535.0 

2H 

835.2 

3H 

1440.5 

4H 

4809.8 

5H 

1503.9 

6H 

945.0 

7H 

760.0 

COMPONENT 

P 

DPP        DHD 

V0L(CU  FT) 

RECUP 

.313      .650 

1.9      EFF=0.91 

TURBINE 

0.7     EFF(P)=0.89 

HEATER 

.019      .430 

31.4 

COMPRESSOR 

10.8      EFF(P)=0.85 

RANKINE  CYCLE:   SATURATION  PRESSURE  =  120.0 
STEAM  FLOW(LBM/S)=  5.69 
FLOW  RATIO  =0.0275 


STATE        TEMP(DEG 

R) 

PRESSURE (PS I) 

IS 

560.0 

0.9 

2S 

575.0 

120.0 

SF 

801.2 

120.0 

3S 

913.0 

120.0 

4S 

560.0 

0.9      QUAL=0.864 

COMPONENT 

VOL(CU  FT) 

BOILER 

507. 

-9 

TURBINE 

6, 

.5 

EFI 

CONDENSER 

134, 

,2 
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V.   APPLICATION 

The  application  of  this  system  to  Navy  ships  has  many 
implications.   One  of  which  is  the  fact  that  neither  Closed 
Brayton  nor  combined  cycles  have  been  previously  used.   However, 
there  is  a  considerable  effort  at  this  time  to  get  COGAS 
(Combined  Gas  Turbine  and  'Steam)  systems  on  board  Navy  ships. 
These  systems  are  all  limited  to  burning  oil  fuels.  The  Closed 
Brayton  combined  cycle  when  equipped  with  a  multi-fuels  heater 
offers  the  additional  flexibility  of  burning  the  cheapest  or 
most  readily  available  fuel. 

Advantages  of  the  Closed  Brayton  combined  cycle  are: 

1)  The  Brayton  cycle  operates  at  temperatures  substan- 
tially below  the  material  limit  than  in  the  open  cycle  gas 
turbines  resulting  in  cheaper,  more  reliable,  and  longer  life 
components. 

2)  Low  pressure  steam  systems  are  well  proven,  and  have 
significantly  less  water  chemistry  problems  than  higher  pressure 
ones . 

3)  As  load  decreases  and  fuel  flow  is  cut  back,  the  con- 
trol valves  (Figure  12)  between  the  compressor  outlet  and  cooler 
inlet  can  be  opened  as  necessary  (keeping  max.  temp  and  speed 
constant)  so  as  to  ensure  that  the  vector  diagrams  of  all  sec- 
tions of  the  turbomachinery  blading  remain  unchanged.   Turbine 
and  compressor  efficiencies  will  remain  essentially  the  same 
producing  high  0  .  at  part  load.   This  requires  use  of  a  CRP 
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Figure   12. 
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(controllable  reversible  pitch)   propeller  to  maintain  constant 
shaft  speed  at  part  loads. 

Disadvantages  of  the  Closed  Brayton  combined  cycle  are: 

1)  Higher  overall  system  weight  and  volume  due  to  the 
additional  heat  exchangers  and  piping  required  for  the  Closed 
Brayton  cycle. 

2)  A  complex,  integrated  control  system  is  required  to 
coordinate  control  of  the  combined  cycle,  fuel  flow,  and  CRP 
at  part  loads. 

3)  Fluidized  bed  combustion  is  not  a  fully  developed 
engineering  concept.   Land  based  prototypes  have  been  built 
and  operated  with  limited  but  increasing  levels  of  success.   A 
shipboard  prototype  is  perhaps  5  years  in  the  future  with  pro- 
duction of  a  marinized  version  at  least  10  years  away. 

4)  If  coal  is  to  be  used  as  the  primary  fuel,  it  should 
be  noted  that  it  has  a  lower  heating  value  per  ton  than  oil. 
The  result  is  a  larger  potentially  more  expensive  (greater 
acquisition  cost)  ship  for  the  same  military  mission. 

5)  Reaction  time  with  coal  fired  fluidized  bed  combustors 
to  load  changes  are  generally  slower  than  that  for  oil  heaters 
necessitating  in  a  separate  fuel  oil  system  for  rapid  increases 
in  load. 
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VI.  CONCLUSIONS  AND  RECOMJV1ENDATIONS 

The  Closed  Brayton  combined  cycle  offers  an  attractive 
alternative  system  for  many  shipboard  applications.  As  with 
the  first  Arab  oil  embargo,  subsequent  cut-offs  of  Arab  oil 
could  result  in  decreased  steaming  days  and  a  lower  overall 
readiness  of  the  U.S.  Navy's  defense  posture. 

If  our  ships  could  have  the  option  of  burning  the  most 
readily  available  fuel,  this  reduction  in  readiness  can  be 
avoided  by  switching  to  coal,  a  fuel  source  in  great  abundance 
in  the  U.S. 

Paramount  in  the  engineering  of  a  successful  multi-fuels 
propulsion  plant  is  the  development  of  a  marinized  fluid  bed 
combustor.   The  design  of  heat  exchangers  and  other  turbo- 
machinery  used  in  the  combined  cycle  are  well  within  the  state 
of  the  art  of  present  technology. 

It  is  thus  that  the  author  concludes  the  following: 

1)  The  design  of  a  modern  coal  burning  combined  cycle  is 
within  the  reach  of  present  technology. 

2)  This  combined  cycle  would  offer  high  efficiency  at 
both  rated  load  and  part  load  due  to  the  unique  controlling 
aspect  of  the  Closed  Brayton  cycle. 

3)  Weight  and  volume  of  a  Closed  Brayton  combined  cycle 
would  be  greater  than  that  of  an  equivalently  rated  COGAS  type 
systemo 

4)  Coal  offers  the  benefit  of  a  readily  available, 
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relatively  cheap  fuel  at  the  expense  of  increased  ship  size 
and  acquisition  cost  due  to  lower  BTU  heat  value  per  ton  when 
compared  to  oil. 

The  author  makes  the  following  recommendations: 

1)  Development  of  a  land  based  test  site  of  a  Closed 
Brayton  combined  cycle  equipped  with  a  fluid  bed  heater  should 
be  considered. 

2)  A  pressurized  fluidized  bed  combustor  suitable  to 
marine  applications  should  be  developed. 

3)  A  shipboard  prototype  system  should  be  integrated 
on  a  U.S.  Navy  ^hip  to  attempt  to  quantitatively  define  the 
shipboard  impact  of  the  system. 


56 


References 

1.  ASME  Steam  Tables,  1967. 

2.  Berman,  I.M.,  "Fluidized  Bed  Combustion  Systems:   Progress 

and  Outlook",  Power  Engineering,  November  1979 • 

3.  "Combined-Cycle  Using  Gas  from  Coal  holds  Promise  for 

Electric  Generation",  Power,  June  1979. 

4.  Critelli,  P.,  Pietsch,  A.  and  Spicer,  N.,  "Closed  Gas 

Turbine  Engines  for  Future  Marine  Propulsion" , 
Marine  Technology,  July  1976. 

5.  Dixon,  S.L.,  Fluid  Mechanics,  Thermodynamics  of  Turbo- 

machinery,  Pergamon  Press,  Oxford,  1966 . 

6.  Foster-?egg,  R.W.,  "A  Coal  Fired  Fluid  Bed  Gas  Turbine 

Cogeneration  System",  Combustion,  June  1979. 

7.  Gruhl,  J.  and  Tearc,  J.D.,  "Empirical  Models  of  Emissions 

and  Energy  Efficiencies  of  Coal -Fired  Fluidized  Bed 
Power  Plants",  M.I.T.  Energy  Laboratory.   Report 
No.  78-015,  Sept.  1977. 

8.  Kavanagh,  Gary  Lee,  "Optimization  of  a  Closed  Brayton 

Cycle  for  Navy  Ship  Conceptual  Design",  Engineer's 
Thesis  in  Ocean  Engineering,  May  1977. 

9.  Keenan,  Joseph  H.  and  Kaye,  Joseph,  Thermodynamic 

Properties  of  Air,  New  York,  John  Wiley  and  Sons,  Inc. 

10.  Lee,  Wing-Hon  Hilary,  "Contribution  to  the  Optimal  Design 

of  Closed  Cycle  Gas  Turbine",  ScD.  Thesis  in  Ocean 
Engineering,  January  1978. 

11.  Muench,  R.K.  and  Knauss,  D.T.  and  Purnell,  J.G.,  "A  Study 

of  Waste-Heat  Boiler  Size  and  Performance  of  a  Con- 
ceptual Marine  COGAS  System",  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center,  February  1980. 

12.  Notes  of  Course  13.21,  "Ship  Power  and  Propulsion"  taught 

by  Dr.  A.  Douglas  Carmichael  at  the  M.I.T.  Department 
of  Ocean  Engineering,  Fall  term  1978. 

13.  Rohsenow,  W.  and  Choi,  H.,  Heat,  Mass .  and  Momentum 

Transfer.  Prentice-Hall,  N.J.,  1961. 


57 


1>.  Roland,  Charles  B.,  "Coal  as  an  Alternative  Fuel  for 
Naval  Vessels",  Engineer's  Thesis  in  Ocean 
Engineering,  May  1978. 

15.  "Subroutine  for  the  Thermodynamic  Properties  of  Steam 

and  Water",  NACSHIPS  0900-009-5010,  Naval  Ship 
Engineering  Center,  Wash.,  D.C.,  1968. 

16.  Vermeir,  Raphael  L. ,  "Combined  Cycles  for  Marine  Propul- 

sions  A  Feasibility  and  Optimization  Analysis" , 
Master's  Thesis  in  Ocean  Engineering,  May  1978. 


58 


APPENDIX  A 

The  computer  program  performs  a  thermodynamic  analysis 
and  component  volume  computation  for  various  specified 
Brayton  cycle  compressor  pressure  ratios. 

The  program  accepts  as  input  the  following  quantities: 

A(P)     -  polytropic  compressor  efficiency  (Brayton 
cycle) 

r, 

/f(P)     -  polytropic  turbine  efficiency  (Brayton 
cycle) 

ER      -  regenerator  effectiveness 

/—^      -   isentropic  pump  efficiency  (Rankine 
cycle) 


(— r-J      -   total  pressure  drop  around  Brayton  cycle 
*  ^  'T  loop 

/IVHB      -   efficiency  of  the  waste  heat  boiler 

T-       -  compressor  inlet  temp.   (Brayton  cycle) 

T^       -   turbine  inlet  temp  (Brayton  cycle) 

T1S      -   steam  condensing  temp. 

AT       -  min.  temp  diff .  between  streams  in  WHB 
(at  P?  and  Tg) 

PSAT     "  m^n'  saturation  pressure  in  the  Rankine  cycle 

Power    -  net  power  out  of  the  combined  cycle 
T.,       -  compressor  inlet  temp.   (heater  cycle) 

T„.      -   turbine  exit  temp,  (heater  cycle) 

•/c^h     -  polytropic  compressor  efficiency  (heater 
cycle) 

/((is)    -   isentropic  turbine  efficiency  (Rankine 
cycle) 
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J^(P)K     _  polytropic  turbine  efficiency  (heater  cycle) 
$  -  equivalence  ratio 

WP       -  working  pressure  in  the  Brayton  cycle. 

Output  is  as  shown  in  Appendix  B. 

The  main  program  performs  the  thermodynamic  analysis  of 
the  closed  Brayton  cycle  and  the  Rankine  cycle.   Functions  per- 
formed by  the  various  subf unctions  are  as  follows: 

1)  Subroutine  HEATER  performs  a  thermodynamic  analysis  of 
the  heater  cycle. 

2)  Subroutine  DHD  computes -^p-  for  the  regenerator,  re- 
cuperator, cooler  and  heater. 

3)  Subroutine  DPP  calculates  ~f~ for  the  regenerator,  re- 
cuperator, cooler,  heater  and  boiler. 

fy)  Subroutine  TURVOL  finds  turbine  volumes. 

5)  Subroutine  COMPVOL  finds  compressor  volumes. 

6)  Subroutine  REGEN  computes  the  regenerator  volume. 

7)  Subroutine  RECUP  computes  the  recuperator  volume. 

8)  Subroutine  COOLER  computes  the  cooler  volume. 

9)  Subroutine  COMBUSTOR  computes  the  heater  volume. 

10)  Subroutine  WHBVOL  computes  the  waste  heat  boiler  volume. 

11)  Subroutine  CONDVOL  computes  the  condenser  volume. 

12)  Subroutine  STMTAB  is  a  computerized  version  of  the  steam 
tables.  Reproduced  from  ref.  15 . 
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DIMENSION  CR(30) ,SATLIQ(3) ,STEAM(3) ,TPP(12) ,WR(11) 
REAL  NCP,NTP,NPIS,NTI$,NWHB,NT,L1,L2,K0 
REAL  NCPH,NTPH,MB,MS,MA,MF,NH,N,IT,KI1,NPP 
REAL  LHV,M,M1,M2,NTU 

COMMON  /STM/PRESS, TEMP, QUAL.SATLIQ, STEAM 
COMMON  /HTR/NH,PRH,DPPH,T2H,T3H,T4H,T5H,T6H,riA,MF 
COMMON  /DIA/DHDREG,DHDCLR,DHDHTR,DHDREC 
COMMON  /PRES/DPPWHB  ,DPPREG  ,DPPCLR  ,DPPREC  ,DPH1  ,DPH2  ,DHDHT 
COMMON  /C0MB/DPPRG,DPPRC,DPPHE,DHDHR,DPHE2 
READ (5, 10)  NCP,NTP,ER,NPIS,NTIS,DP,NWHB,T1,T4,T1S,DT,PMIN,K 
FORMAT  (7F4.3,5F7.1,I2) 
READ (5, 15)  PWR,T1H,T7H,NCPH,NTPH,PHI,WP 
FORMAT  (F8.1,2F6.1,3F5.3,F6.1) 
READ(5,20)  (CR(I),I=1,K) 
FORMAT  (<K>F6.2) 
R=  .068544 
DO  200  1=1, K 
PR=CR( I ) 

Z=CP1(T1)+(R/NCP)*L0G(PR) 
T2=EXP(Z/.2475) 
CALL  ITER  (T2,Z) 
TR=1/(1-DP)/PR 
Z=CP1(T4)+R*NTP*L0G(TR) 
T5=EXP(Z/.2475) 
CALL  ITER  (T5,Z) 
IF(T2  .GT.  T5)  GO  TO  30 
T6=T5-ER*(T5-T2) 
T3=T2+ER*(T5-T2) 
GO  TO  40 
T6=T5 
T3=T2 
L1=T1 
L2=1122.0 

IF(L2  .GT.  T6)  GO  TO  50 
GO  TO  60 
50     L2=T6-(T6-T1)/10. 
60     TPP(1)=L2 
WR1=0. 

C=(L2-L1)/10. 
DO  70  J=l,ll 
TSF=TPP(J)-DT 
T3S=T6-DT 
T=TSF-460 

CALL  STMTAB  (2,T,0.) 
PSF=PRESS 
HSF=SATLIQ(2) 
T=T3S-460 

CALL  STMTAB  (7,PSF,T) 
H3S=STEAM(2) 
S3S=STEAM(3) 
T=TlS-460 
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CALL  STMTAB    (2,T,0.) 

P4S=PRESS 

V1S=SATLIQ(1) 

H1S=SATLIQ(2) 

S1S=SATLIQ(3) 

CALL  STMTAB   (9,P4S,S3S) 

H4SIS=STEAM(2) 

H4S=H3S-NTIS*(H3S-H4SIS) 

CALL  STMTAB    (8,P4S,H4S) 

X=QUAL 

V4S=STEAM(1) 

H2SIS=H1S+V1S*(PSF-P4S) 

H2S=H1S+(H2SIS-H1S)/NPIS 

T2S=T1S+15. 

B=NWHB*(CP2(T6)-CP2(TPP(J)))/(H3S-HSF) 

Z=ACP(TPP(J))*TPP(J)-B*(HSF-H2S)/NWHB 

T7=TPP(J) 

CALL  ITER2    (T7,Z) 

WR(J)=B*(H3S+H1S-H4S-H2S) 

IF(J    .EQ.    1)   GO  TO  68 

IFCT7    .LT.   Tl)  GO  TO  75     ' 

IF(PSF   .LT.   PMIN)   GO  TO  75 

IF(WR(J)    .GT.   WR1)  GO  TO  68 

IF(X    .GT.    .87)   GO  TO  75 
68  WR1=WR(J) 

TPP(J  +  1)=TPP(J)-C 
70      CONTINUE 
GO  TO  76 

75  IF(ABS(TPP(J)-TPP(J-D)  .LT.  .1)  GO  TO  76 
L1=TPP(J) 

L2=TPP(J-1) 
GO  TO  60 

76  WB=CP2(T4)+CP2(T1)-CP2(T5)-CP2(T2) 
QIN=CP2(T4)-CP2(T3) 
NT=(WB+WR(J))/QIN 
MB=.707*PWR/(WB+WR(J)) 

MS=B*MB 

QB=MB*QIN 

CALL  HEATER  (QB  ,T1H,T7H,NCPH,NTPH,PHI,T3 ,T4) 

CALL  DHD(PR) 

DHDHT=DHDHTR 

K=l 

77  CALL  DPP(T6,T7,PR,DP,ER,T2,T3,DPPH,T4H,T5H,T4,PRH,WP,MB,MA 
*  ,MF,K) 

G=l  33 

N=G/(G-(G-1.)*NTP) 

P=WP*(PR-DPPREG) 

V=1./ARH0(T4,P) 

VMNTURB=0.0 

CALL  TURBVOL  (N,MB,V,TR,T4  .VMNTURB) 

G=1.39 
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N=G/(G-(G-1.)*NTPH) 

P=14.7*(PRH-DPPREC) 

V=1./ARH0(T6H,P) 

M=MA+MF 

VHTRTURB=0.0 

TRH=1./(1.-DPPH)/PRH 

CALL  TURBVOL   (N,M,V,TRH,T6H,VHTRTURB) 

T=T1S/T3S 

P=P4S/PSF 

N=l./(l.-LOG(T)/LOG(P)) 

VSTMTURB=0.0 

CALL  TURBVOL   (N,MS,V4S,P  ,T3S,VSTMTURB) 

G=1.4 

N=G/(G-(G-1.)/NCP) 

V=l./ARHO(Tl,WP) 

VMNCOMP=0.0 

CALL  COMPVOL    (N,MB  ,V  ,PR,VMNCOMP) 

V=1./ARH0(535. ,14.7) 

VHTRCOMP=0.0 

CALL  COMPVOL    (N,MA,V,PRH,VHTRCOMP ) 

D=.125 

VREGEN=0.0 

CALL  REGEN   (T2,T3,T5,T6  ,WP  ,D,MB,PR,DPPREG,DHDREG,VREGEN) 

VRECUPO.O 

CALL  RECUP   (T2H,T3H,T5H,T6H,PRH,DPPH,D,MA,MF,DPPREC,VRECUP) 

D=.25 

VCLR=0.0 

CALL  COOLER  (T7  ,T1  ,WP  ,DSNTP,MB ,DHDCLR,DPPCLR,VCLR) 

DPPRG=DPPREG 

DPPRC=DPPREC 

DPPHE=DPPH 

DHDHR=DHDHTR 

DPHE2=DPH2 

VHTR=0.0 

CALL  COMBUSTOR  (T3,T4,T4H,T5H  ,WP  ,PR,D,MB,PRH,MA,MF,VHTR) 

VWHB=0.0 

CALL  WHBVOL   (MB  ,MS,T6  ,T7  ,NWHB,H2S,H3S,DPPWHB  ,VWHB) 

VC0ND=0.0 

CALL  CONDVOL   (MS,X,VCOND) 

VOLT=VMNTURB+VHTRTURB+VSTMTURB+VMNCOMP+VHTRCOMP+VREGEN+VRECUP 
*     +VCLR+VHTR+VWHB+VCOND 

IF(K   .EQ.    1)  GO  TO  78 

IF(VOLT   .GT.   VOL1)   GO  TO  81 

IFUDPPCLR-.002)    .LT.   O.O)  GO  TO  81 
78  K=K+1 

V0L1=V0LT 

V1=VCLR 

V2=VWHB 

GO  TO  77 
81  V0LT=V0L1 

VCLR=V1 
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VWHB=V2 

WRITE(6,82)   PR, WP, DP, MB 

82  FORMAT (  '1','BRAYTON  CYCLE:   PRESSURE  RATI 0= ' ,F5.2/1X, 'WORKING 

*  PRESSURE^, F5.1/1X,  'TOTAL  PRESSURE  DROP= '  ,F3.2/1X, 'MASS  FLOW 

*  (lbm/s)=',F6.2//10X, 'STATE     TEMP(deg  R) ' ) 

WRITE(6,83)   Tl 

83  FORMAT (12X,'l' ,6X,F6.1) 
WRITE(6,84)   T2 

84  F0RMAT(12X,'2',6X,F6.1) 
WRITE(6,85)   T3 

85  F0RMAT(12X,'3' ,6X,F6.1) 
WRITE(6,86)   T4 

86  FORMAT (12X, '4' ,6X,F6.1) 
WRITE(6,87)   T5 

87  F0RMAT(12X,'5' ,6X,F6.1) 
WRITE(6,88)   T6 

88  F0RMAT(12X,'6',6X,F6.1) 
WRITE(6,89)   TPP(J) 

89  F0RMAT(12X,'PP' ,5X,F6.1) 
WRITE(6,90)   T7 

90  F0RMAT(12X,'7',6X,F6.1//5X, 'COMPONENT     DPP     DHD     VOL(cuft)') 
WRITE (6, 91)  DPPREG,DHDREG,VREGEN,ER 

91  F0RMAT(7X,'REGEN' ,4X,F4.3,1X,F4.3,1X,F7 .1 ,2X, 'EFF=' ,F4 .2) 
WRITE(6,92)   VMNTURB,NTP 

92  F0RMAT(6X, 'TURBINE' ,15X,F5 .1,2X, 'EFF(P)= 'F4.2) 
WRITE(6,93)   DPPCLR,DHDCLR,VCLR 

93  F0RMAT(6X, 'COOLER' ,4X,F4.3,1X,F4.3,3X,F5 .1) 
WRITE (6, 94)   VMNCOMP,NCP 

94  FORMATUX,  'COMPRESSOR'  ,14X,F5  .1 ,2X,  'EFF(P)=  '  ,F4.2) 
WRITE(6,95)  DPH1 

95  F0RMAT(6X,'HEATER' ,4X,F5.4) 
WRITE(6,96)  DPPWHB,NWHB 

96  F0RMAT(6X, 'BOILER' ,4X,F4.3,15X, 'EFF= ' ,F4.2/) 
WRITE(6,100)   PRH,DPPH,MA,MF 

100  F0RMAT(1X, 'HEATER  CYCLE:   PRESSURE  RATIO= ' ,F5 .2/1X, 'TOTAL 

*  PRESSURE  DR0P=',F3.2/1X,'AIR  FLOW(lbm/s)=' ,F6.2/1X, 'FUEL 

*  FL0W(lbm/s)  =  ',F6.2//10X,'STATE     TEMPtdeg  R) ' ) 

WRITE(6,101)   T1H 

101  F0RMAT(12X,'lh' ,5X,F6.1) 
WRITE(6,102)   T2H 

102  F0RMAT(12X,'2h'  ,5X,F6.1) 
WRITE(6,103)   T3H 

103  F0RMAT(12X,'3h' ,5X,F6.1) 
WRITE(6,104)   T4H 

104  F0RMAT(12X,'4h' ,5X,F6.1) 
WRITE(6,105)   T5H 

105  F0RMAT(12X,'5h' ,5X,F6.1) 
WRITE(6,106)   T6H 

106  F0RMAT(12X,'6h' ,5X,F6.1) 
WRITE(6,107)   T7H 

107  F0RMAT(12X,'7h' ,5X,F6 .1//5X, 'COMPONENT     DPP     DHD     VOL(cu   ft) ' 


6^ 


*  ) 

ERH= ( T3H-T2H ) / (T5H-T2H  ) 

WRITE (6 ,110)  DPPREC ,DHDREC ,VRECUP,ERH 

110  FORMAT! 7 X, 'RECUP' ,4X,F4.3,1X,F4.3,1X,F7.1,2X, 'EFF='  ,F4.2) 
WRITE(6,111)  VHTRTURB,NTPH 

111  F0RMAT(6X,' TURBINE', 15X,F5 .1 ,2X, 'EFF(P)= ' ,F4. 2) 
WRITE(6,112)  DPH2,DHDHTR,VHTR 

112  F0RMAT(6X, 'HEATER' ,4X,F4.3,1X,F4.3,3X,F5.1) 
WRITE(6,113)  VHTRC0MP,NCPH 

113  F0RMAT(4X, 'COMPRESSOR' ,14X,F5.1,2X, 'EFF(P )  =  ' ,F4.2/ ) 
WRITE(6,120)  PSF,MS,B 

120  F0RMAT(1X,'RANKINE  CYCLE:  SATURATION  PRESSURE= ' ,F6.1/1X, 

*  'STEAM  FL0W(lbm/s)=',F6.2/lX,'FL0W  RATI 0=' ,F6. 4/ /10X, 'STATE 

*  TEMPI deg  R)  PRESSURE(psi) ' ) 

WRITE(6,121)  T1S,P4S 

121  F0RMAT(12X,'ls' ,5X,F6.1 ,8X,F6 .1 ) 
WRITE(6,122)  T2S,PSF 

122  F0RMAT(12X,'2s',5X,F6.1,8X,F6.1) 
WRITE(6,123)  TSF,PSF 

123  F0RMAT(12X, 'SF' ,5X,F6 .1,8X,F6 .1) 
WRITE(6,124)  T3S,PSF 

124  F0RMAT(12X,'3s',5X,F6.1,8X,F6.1) 
WRITE(6,125)  T1S,P4S,X 

125  F0RMAT(12X,'4s'  ,5X,F6  .1 ,8X,F6  .1 ,2X, 'QUAL=  '  ,F5  .3/) 
WRITE(6,128)  VWHB 

128  F0RMAT(5X, 'COMPONENT     VOL(cu  ft)  76X, 'BOILER' ,5X,F6  .1) 
WRITE(6,129)   VSTMTURB,NTIS 

129  F0RMAT(6X, 'TURBINE', 5X,F5.1,2X,'EFF(IS)=',F4. 2) 
WRITE (6, 130)   VCOND 

130  F0RMAT(5X, 'CONDENSER' ,4X,F5.1/) 
WRITE(6,135)   NT 

135  FORMATCl' , 'SUMMARY: THERMAL  EFF=',F5.3) 
WRITE(6,136)   NH,PHI 

136  F0RMAT(9X, 'HEATER  EFF=',F5.3,'      (PHI= ' ,F4.2, ' ) ' ) 
SFC=MF*3600./PWR 

WRITE(6,137)   SFC,VOLT,PWR 

137  F0RMAT(9X,'SFC(lbm/HP-hr)=',F5.3/9X,'V0L  TOTAL(cu  ft)=',F6.1/ 
*9X,'P0WER(HP)  =  ',F7.1) 

RANKINE=WR(J)/WB 
WRITE(6,138)  RANKINE 

138  FORMAT (9X, 'WORK  FRAC  BY  STM  CYCLE= ' ,F5 .3) 
200     CONTINUE 

-  STOP 
END 

FUNCTION  CPKT) 

CP1=.2475*L0G(T)-(3.759E-5*T)+(2.553E-8*T**2)-(4.103E-12*T**3) 
RETURN 
END 

FUNCTION  CP2(T) 

CP2=.2475n-(1.879E-5*T**2)  +  (1.702E-8n**3)-(3.078E-12*T**4) 
RETURN 
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END 

FUNCTION  AMU(T) 
TK=5.*T/9 . 

AMU=TK**l.*5*l.E-6/(TK+110.4) 
RETURN 
END 

FUNCTION  ARHO(T,P) 
R=53.35 

ARH0=P*144./R/T 
RETURN 
END 

FUNCTION  AK(T,P) 
A=.6325E-5 
B=245.4 
C=12. 

TK=5.*T/9 . 

AKO(A*TK**.5)/(1.+B*10.**(-1.*C/TK)/TK)*418.6 
AK=(  AK0+U.464E- 5*  (16. 02*ARH0(T,P))**1. 23))*. 578/3600. 
RETURN 
END 

FUNCTION  APR(T,P) 
CP=ACP(T) 
XMU=AMU(T) 
XK=AK(T,P) 
APR=CP*XMU/XK 
RETURN 
END 

FUNCTION  ASIGMA(T,P) 

ASIGMA=17.214*(6228.*AK(T,P))**(-1.41)*(1.46*AMU(T))**1.23* 
*(1./APR(T,P)**.564)*(1./(16.02*ARH0(T,P))**.41) 
RETURN 
END 

FUNCTION  ACP(T) 

ACP=.  2475-3. 759E-5*T+5.106E-8*T**2-1.231E-U*T**3 
RETURN 
END 

SUBROUTINE  ITER   (T,X) 
10       T1=T 

T=EXP((X+.2475*L0G(T)-CPl(T))/.2475) 

IF(ABS(T-T1)    .LT.    l.E-2)   GO  TO  20 

GO  TO  10 
20       RETURN 

END 

SUBROUTINE  ITER2(T,X) 
10       T1=T 

T=X/ACP(T) 

IF(ABS(T-T1)    .LT.   l.E-2)   GO  TO  20 

GO  TO  10 
20       RETURN 

END 
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SUBROUTINE  HEATER   (QB ,T1H,T7H,NCPH,NTPH,PHI ,T3,T4) 

REAL  LHV,NH,MF,MA,NCPH,NTPH,N 

COMMON   /HTR/NH,PRH,DPPH,T2H,T3H,T4H,T5H,T6H,MA,MF 

LHV=12900. 

FA=.0922 

H7H=ACP(T7H)*T7H 

H1H=ACP(T1H)*T1H 

NH«1-((1+FA*PHI)*H7H-H1H)/(FA*PHI*LHV) 

MF=QB/(NH*LHV) 

MA=MF/(FA*PHI) 

T2H=0.0 

T6H=0.0 
10       CPT=ACP((T6H+T7H)/2.) 

CPC=ACP((TlH+T2H)/2.) 

N= (1+FA*PHI )*CPT*NTPH* ( CPT*T7H-CPC*T1H ) 

D=CPC*T1H*(CPC/NCPH-CPT*NTPH) 

PRH=(N/D)**(NCPH/CPC) 

TH2=T1H*PRH**(CPC/NCPH) 

N=(1+FA*PHI)*CPT*T7H+CPC*T1H*(PRH**(CPC/NCPH)-1) 

D=(1+FA*PHI)*CPT*T7H*PRH**(CPT/NTPH) 

DPPH=1-(N/D)**(1./(CPT*NTPH)) 

TRH=(1-DPPH)*PRH 

TH6=T7H*(TRH)**(CPT*NTPH) 

IF(ABS(T2H-TH2)    .LT.    .1)   GO  TO  20 

IF(ABS(T6H-TH6)    .LT.    .1)   GO  TO  20 

T2H=TH2 

T6H=TH6 

GO  TO  10 
20       T2H=TH2 

T6H=TH6 

T5H=T3+100. 

T4H=T5H+QB/(ACP(T5H/2+1300. )*(MA+MF) ) 

T3H=T2H+ ( 1+MF/MA)* ( T5H-T6H ) 

RETURN 

END 

SUBROUTINE  DHD(PR) 

COMMON   /DIA/DHDREG,DHDCLR,DHDHTR,DHDREC 

DHDREG=.557+.0933*PR 

IF(DHDREG   .GT.    .95)  DHDREG=.95 

DHDCLR=.076+.092*PR 

DHDHTR=.214+.048*PR 

DHDREC=.65 

RETURN 

END 

SUBROUTINE  DPP(T6  ,T7  ,PR,DP,ER,T2,T3,DPPH,TH4,TH5,T4,PRH,P  ,M,MA 
*,MF,K) 
REAL  M,MA,MF,M1,M2 
COMMON  /PRES/DPPWHB,DPPREG,DPPCLR,DPPREC,DPH1,DPH2,DHDHT 

67 


DTG=T6-T7 

DPPWHB=FL0AT(K)*.Q02 

DPPREG=( .135-.0284*PR)*DP/.l 

IF(DPPREG    .LT.    .01)  DPPREG=.01 

IF(T2    .EQ.   T3)   DPPREG=0.0 

DPH2=(.0108*PR-.01)*DP/.l 

T1=(T3+T4)*.5 

TT=(TH4+TH5)*.5 

P1=P*(PR-DP) 

P2=14.7*(PRH-DPPH) 

M1=M 

M2=MA+MF 

XMU1=AMU(T1) 

XMU2=AMU(TT) 

XRH01=ARH0(T1,P1) 

XRH02=ARH0(TT,P2) 

B=(M1/M2)**(-1.75)*(XMU1/XMU2)**(-.25)*(XRH01/XRH02)*(P1/P2) 

DPH1=DPH2*DHDHT**3/B 

DPPCLR=DP-DPPWHB-DPPREG-DPH1 

DPPREC=DPPH-DPH2-.l 

RETURN 

END 

SUBROUTINE  TURBVOL   (G,M,V  ,TR,T,V0L) 
REAL  M,N,KI1,IT 
SF=4.0 
UT=850. 
PHIT=0.6 
VA=0.9 
CL=2.0 
PI=3. 14159 
RH0M=15.12 
A=9.63E6 
B=-95.14 
ALPHA=1.2 
Al=2.5 
Xl=.307 
KIl=3.39E-3 
THETA=0.8 

RT=SQRT(M*V/(PHIT*UT*PI*(1.-VA**2))) 
VB*SQRT(1.-TR**(1./G)*(1.-VA**2)) 
C1=0.2*(A1+1.)*1./V*(1.-VA**2)*UT**2*RT**2*X1/KU 
C2=2.*(A-B*T)/(SF*ALPHA) 
C3=2.*B*T*(1.-VA**2)**(G-1.)/(SF*ALPHA) 
C4«-l.*THETA*RH0M*UT**2/2. 
H=(VA-VB)/20. 
SUM=0.0 
X1=0.0 
Y1=VA 

DO  100  J=l,21 
N=C1*((1.-Y1)/(1.+Y1))*PHIT**2+Y1**2 
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D=C2+C3*(1.-Y1**2)**(1.-G)+C4*(1.-Y1**2) 

DN=C1*(2.*Y1*(1.-Y1)/(1.+Y1)-(PHIT**2+Y1**2)/(1.+Y1)-(1.+Y1) 
***(-2)*(l.-Yl)*(PHIT**2+Yl**2)) 

DD=2.*Y1*(1.-Y1**2)**(-1.*G)*(G-1.)*C3-2.*Y1*C4 
X2=.5*(1./(N/D)**.5)*((DN*D-DD*N)/D**2) 
IF(X1  .EQ.  0.0)  GO  TO  50 
SUM=SUM+(Xl+X2)/2. 
50     X1=X2 

Y1=Y1-H 
100   CONTINUE 
IT=H*SUM 

V0L=2.*PI*RT**2*CL*IT 
RETURN 
END 

SUBROUTINE  COMPVOL   (G,M,V,PR,VOL) 

REAL  M,N,KI1,IT 

SF=4.0 

UT=850. 

PHIT=0.5 

VA=0.6 

CL=2.0 

PI=3. 14159 

RH0M=15.12 

ALPHA=1.2 

Al=12.0 

Xl=.05 

KIl=4.93E-5 

THETA=0.8 

SIGMAY=12.54E6 

RT=SQRT(M*V/(PHIT*UT*PI*(1.-VA**2))) 

VB*SQRT(1.-PR**(-1./G)*(1.-VA**2)) 

C1=.2*(A1+1.)*1./Y*(1.-VA**2)*UT**2*RT**2*PHIT**2*X1AI1 

C2=-l*THETA*RH0M*UT**2/2. 

C3=2.*SIGMAY/(SF*ALPHA) 

H=(VB-VA)/20. 

SUM=0.0 

X1=0.0 

Y1=VA 

DO  100  J=l,21 

N=C1*((1.-Y1)/(1.+Y1))*(PHIT**2  +  .49*Y1**2)/(PHIT**2  +  .7225*Y1**2 
*) 

D=C3+C2*(1.-Y1**2) 

DN=C1*(-1.*(PHIT**2  +  .49*Y1**2)/(PHIT**2  +  .7225*Y1**2)*(1.+Y1)-(1 
*.-Yl)*(l.+Yl)**(-2)*(PHIT**2+.49*Yl**2)/(PHIT**2+.7225*Yl**2)  + 
*.98*Y1*(1.-Y1)/(1.+Y1)/(PHIT**2+.7225*Y1**2)-1.445*Y1*(1.-Y1)/ 
*(l.+Yl)*(PHIT**2+.49*Yl**2)*(PHIT**2+.7225*Yl**2)**(-2)) 

DD=-2.*C2*Y1 

X2=.5*(1./(N/D)**.5)*((DN*D-DD*N)/D**2) 

IF(X1    .EQ.   0.0)   GO  TO  50 

SUM=SUM+(Xl+X2)/2. 
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50   X1=X2 

Y1=Y1+H 
100  CONTINUE 

IT=-H*SUM 

V0L=2.*PI*RT**2*CL*IT 

RETURN 

END 

SUBROUTINE  REGEN  (T2,T3,T5  ,T6  ,P,D,MB,PR,DPPREG,DHDREG,VOL) 
REAL  MB 
T=T5/2+T6/2 
P1=P*PR 

SIGMA=ASIGMA(T,P1) 
IF(T2  .EQ.  T3)  GO  TO  99 
DTLM=T3-T2 

GAMMA=1./(6803.*P1)**.41*(.0254*D)**1.18/(.56*DTLM)**1.41 
DH1=CP2(T5)-CP2(T6) 
PI=3. 14159 

AREG=PI/4.*SIGMA*GAMMAMMB*.454)*(DH1*2326J**1.41*(1.+DHDREG) 
***2.41*(1.+PR**(+2)*DHDREG**(-3))**.41*35.3 
V0L=AREG*DPPREG**(-.41) 
99   RETURN 
END 

SUBROUTINE  RECUP  (T2H  ,T3H,T5H  ,T6H,PRH,DPPH,D,MA,MF,DPPREC,V0L) 
REAL  MA,MF 
T=T5H/2.+T6H/2. 
PC=1./(1.-DPPH) 
P=14.7*PRH 
SIGMA=ASIGMA(T,P) 
DTLM=T3H-T2H 

GAMMA=(6803.*P)**(-.41)*(.0254*D)**1.18*{.56*DTLM)**(-1.41) 
DH1=CP2(T5H)-CP2(T6H) 
PI=3. 14159 
DHDREC=.65 
CM=MA+MF 

AREC=PI/4.*SIGMA*GAMMA*(CM*.454)^DH1*2326.)**1.41*(1.+DHDREC) 
***2.41*{1.+0HDREC**(-3)*PC**2)**.41*35.3 
V0L=AREC*DPPREC**(-.41) 
RETURN 
END 

SUBROUTINE  COOLER  (T7  ,T1,P,D,NTP,MB,DHDCLR,DPPCLR,V0L) 

REAL  NPP,NTP,MB 

IF(T7    .LT.  Tl)   GO  TO  99 

TA=T7/2.+Tl/2. 

SIGMA=ASIGMA(TA,P) 

TC1=530. 

TC2=540. 

DTLM=(T1-TC1-T7+TC2)/L0G((T1-TC1)/(T7-TC2)) 

GAMMA=(6803.*P)**(-.41)*( .0254*D)**1.18*( .56*DTLM)**(-1.41 ) 
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DH1=CP2(T7)-CP2(T1) 
PI=3. 14159 
TW=TCl/2.+TC2/2. 
RK=AK(TA,P)/l.E-4 
RMU=AMU(TA)/6.4E-4 
RM=(T7-T1)*ACP(TA)/(1.*10.) 
RPR=APR(TA,P)/6.4 
NPP=.80 

RRH0=ARH0(TA,P)/64.0 
RSV=ARH0(T7,P)/64.0 
ALPHA=RK*RMU**(-.8)*RM**.8*RPR**.4 
CB=RM**(-2.75)*RMU**(-.25)*RRH0*RSV/NTP/NPP 

ACLR=PI/4.*SIGMA*GAMMA*(MB*.454)*(DH1*2326.)**1.41*(1.+DHDCLR) 
**(1.+ALPHA*DHDCLR)**1.41*(1.+CB*DHDCLR**(-3))**.41*35.3 
V0L=ACLR*DPPCLR**(-.41) 
99   RETURN 
END 

SUBROUTINE  COMBUSTOR  (T3,T4,T4H,T5H,WP  ,PR,D,MB,PRH,MA,MF  ,VOL) 
REAL  MB,MA,MF,M 

COMMON  /C0MB/DPPRG,DPPRC,DPPHE,DHDHR,DPHE2 
T=T3/2+T4/2 
P=WP*(PR-DPPRG) 
SIGMA=ASIGMA(T,P) 

DTLM=(T5H-T3-T4H+T4)/L0G((T5H-T3)/(T4H-T4)) 
GAMMA=(6803.*P)**(-.41)*(.0254*D)**1.18*(.56*DTLM)**(-1.41) 
DH1=CP2(T4)-CP2(T3) 
PI=3. 14159 
TH=T4H/2+T5H/2 
TRH=(1.-DPPHE)*PRH 
PH=14.7*(PRH-DPPRC) 
RK=AK(T,P)/AK(TH,PH) 
RMU=AMU(T)/AMU(TH) 
M=MA+MF 
RM=MB/M 

RPR=APR(T,P)/APR(TH,PH) 
ALPHA=RK*RMU**(-.8)*RM**.8*RPR**.4 

AHTR=PI/4.*SIGMA*GAMMA*(MB*.454)*(DH1*2326.)**1.41*(1.-K)HDHR)* 
*     (l.+ALPHA*DHDHR)**1.41*DHDHR**(-1.23)*35.3 
V0L=AHTR*DPHE2**(-.41) 
RETURN 
END 

SUBROUTINE  WHBVOL  (MB,MS,T6  ,T7  ,NWHB,H2S  ,H3S  ,DPPWHB,YOL) 

REAL  MB,MS,NWHB,NTU 

T=T6/2+T7/2 

CP=ACP(T) 

CPS=2.0 

CR=MS*CPS/(MB*CP) 

IF(CR  .GT.  1.)  GO  TO  10 

CMIN=MS*CPS 

71 


GO  TO  20 
10   CR=1./CR 

CMIN=CP*MB 
20   PS=4.0 

ER*(NWHB*CR-1.)/(NWHB-1.) 

NTU=(-PS)*L0G(1.+1./CR*L0G(1.-CR*(1.-ER**(1./PS))/(CR-ER**(1./ 
*PS)))) 

DH=H3S-H2S 

VOL=.00045*CMIN*NTU*MS*DH*DPPWHB**(-.41) 

RETURN 

END 

SUBROUTINE  CONDVOL  (MS,X,VOL) 

REAL  MS 

V0L=MS*X*27.3 

RETURN 

END 
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SUBROUTINE  STMTAB 
DIMENSION  SATLIQ( 
COMMON  /STM/PRESS 


DATA  (COEFFT(I) 
.69573391E+03 

-.19513682E+00 

-.50625849E+00 
.55592942E+03 
.26548914E-04 
.30706919E+00 
.45000000E+03 
DATA  (COEFFT(I) 

-.32893178E-14 
.81536841E-02 

-.64272696E-21 

-.64251564E-10 

.33953897E+01 

.92756543E-17 

.39936598E-08 

DATA  (COEFFT(I) 

-.14722587E-21 
.27455505E-12 
.98512754E+00 
.73857180E-22 

-.70679233E-14 
.30734203E-04 

-.33608032E-23 
DATA  (COEFFT(I) 
.15587770E-07 
.55000000E+03 
.44979672E-12 
.12399614E-02 
.12945371E-21 
.17548001E-12 
.10150974E+01 
DATA  (COEFFT(I) 
.69217855E-16 
.59598410E-08 
.25000000E+03 
.49239634E-17 

-.85546496E-06 
.19000000E+01 

-.18796109E+03 
DATA  (COEFFT(I) 
.11199490E+04 
.62379937E+05 
.35122466E+04 
.30000000E+01 
.52122736E+03 
.12340954E-04 
.18000000E+01 


(NI.XI.TI) 

3),STEAM(3)5C0EFFT(237) 

, TEMP, QUAL,SATLIQ, STEAM 

=1,30)  /    .22082221E+04,    .65000000E+03 

.60464261E+01,    .70000000E+01,    .18552680E+00 


.52279952E+00 

.72595720E+00 

.20000000E+01 

.91349306E-04 

.36241343E+01 

.55000000E+03 

=  31,  60) 

.11734466E-11 

.69519265E+01 

.12704366E-18 

.44797322E-07 

.43010158E+03 

.95470656E-15 

.99567313E-06 

=  61,  90) 

.66553573E-19 

.18335302E-10 

.25000000E+03 

.33225506E-19 

.90114201E-11 

.21762691E-01 

.85819792E-20 

=  91,120) 

.81264440E-05 

.94743147E-18 

.51337380E-10 

.12798551E+01 

.65873139E-19 

.23202424E-09 

.21848036E+03 

=121,150) 

.18345822E-14 

.45575915E-06 

.50084539E-23 

.16250351E-13 

.14285709E-02 

.11085077E+05 

.15699797E+03 

=151,180) 

.15000000E+01 

.29474498E+05 

.12340833E+04 

.13000000E+04 

.11682947E-01 

.13555136E-03 

.23822091E+02 


.91148958E-01 

.14585694E+02 

-.74133440E-06 

.15679922E-02 

.44960616E+02 

-.98736888E-22 

/  .30041189E-17 

.11244476E-07 

.25000000E+03 

-.96976163E-16 

-.28104756E-04 

.54927377E+03 

-.29004677E-12 

-.59157992E-03 

/    .55000000E+03 

.47725605E-17 

-.31160333E-06 

.47000000E+03 

.30075575E-18 

.58149642E-09 

.25000000E+03 

.25729365E-18 

/    .60506309E-13 

.17002639E-01 

.24971478E-16 

.49178986E-08 

.54927377E+03 

-.13817574E-17 

.36145967E-06 

.45000000E+03 

/-.21134438E-20 

-.66417458E-12 

.12325099E-02 

-.15659070E-21 

-.74068320E-12 

.36753405E+00 

-.15994587E+05 

-.20213342E+03 

/  .15040455E+03 

-.14863585E+06 

-.18907586E+05 

-.68928827E+02 

.70848990E+01 

-.12001588E+00 

-.15099319E-03 

.18776649E+03 


.61843682E-01 

.12404319E+03 

.73324823E-06 

.25032530E-01 

.17926732E+03 

.19301230E-19/ 

-.42113669E-15 

-.78419653E-05 

.20575548E-23 

.83370000E-13 

.17490541E-01 

-.18444779E-19 

.35353292E-10 

.78143274E+00/ 

.21848036E+03 

-.53349163E-14 

-.11956354E-03 

.55000000E+03 

-.73898651E-15 

.86712981E-07 

-.78459129E-25 

-.31591904E-15/ 

.13684599E-10 

.45000000E+03 

-.38302588E-14 

.48326895E-05 

.25000000E+03 

.62520106E-14 

.12506461E-03 

.55000000E+03/ 

.15327038E-18 

.14589563E-10 

.74969797E+00 

-.46099119E-18 

.10479899E-08 

.17000000E+01 

.91407688E+04 

-.48650176E+02/ 

-.19050547E+03 

-.13414819E+06 

-.57858876E+03 

.12036836E+04 

.16507788E+02 

.20411991E+01 

.30000000E+01 

.33080504E+03/ 
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DATA  (COEFFT(I), 1=181, 210)          /- 

.32141408E+01, 

.42449101E+03 

*,    .16400835E+04,    .18746657E+02, 

.63978878E+02, 

.20719841E+04 

*,    .13000000E+04,    .18000000E+01,- 

•.64369800E-05, 

.96418803E-02 

*,    .38879248E+01,    .11034060E-05, 

.47225001E-03,- 

.94113600E+01 

*,    .69416973E-05,-.60892120E-02, 

.10951922E+01, 

.13000000E+04 

*,    .18000000E+01,-.33101208E-03, 

.21661916E+01, 

.53500291E+03 

*,-.77236618E-03,    .36753650E+00,- 

.25450121E+03, 

.25249707E-02 

*,-.16427106E+01,    .61477766E+03, 

0.0, 

0.0           / 

DATA  (COEFFT(I), 1=211, 237)          / 

0.0,   0.0, 

0.0,   0.0,   0.0 

*,  0.0,0.0,0.0,0.0,0.0,0.0,0.0,0. 

0,0.0,0.0,0.0, 

.90000000E+03 

*,    .14000000E+01,-.15478259E-04, 

.19351556E-01, 

.24698909E+01 

*,    .23329644E-04,    .11493791E-02,- 

.13139208E+02, 

.58860219E-05 

*,-.12592899E-01,    .55101160E+01/ 

N=NI 

X=XI 

Y=YI 

I=N-3 

LINK=0 

PRESS=X 

DELTA1=0.0 

ERR0R1=0.0 

ERR0R2=0.0 

IF  (N-6)   1000,1000,1030 

1000  IF   (I)   1020,1020,1010 

1010  LINK=1 

N=I 

1020  IF   (X-COEFFT(N))   1040,1040,1620 

1030  N=I 

1040  GO  TO   (1130,1060,1190,1230,1080,1090,1180,1330) 

,N 

1050  X=Z 

N=2 

1060  TEMP=X 

GO  TO  1190 

1070  PRESS=EXP(Z) 

GO  TO  1150 

1080  11=2 

JJ=3 

T0LER=.005 

GO  TO  1100 

1090  11=3 

JJ=2 

TOLER=  .000005 

1100  H=Y 

1110  LINK=2 

1120  N=l 

1130  X=ALOG(PRESS) 

PLOG=X 

GO  TO  1190 

1140  TEM>=Z 

1150  X=TEMP 

J=l 
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N=4 

GO  TO  1190 
1160  SATLIQ(J)=Z 

J=J+1 

N=N+1 

IF  (J-3)  1190,1190,1170 
1170  IF  (LINK)  1640,1640,1300 
1180  H=X 

S=Y 

x=s 

1190  I=23*N-19 

K=l 

IF  (X-COEFFT(I))  1200,1200,1210 
1200  K=12 
1210  I=I+K 

ARG=X-COEFFT(I) 

Z=0.0 

DO  1220  K=l,10 

1=1+1 
1220  Z=Z*ARG+COEFFT(I) 

GO  TO  (1140,1070,1050,1160,1160,1160,1410) ,N 
1230  TEMP=Y 

GO  TO  1300 
1240  ERROR=H-STEAM(II) 

IF   (ABS(ERROR)-TOLER)   1640,1250,1250 
1250  IF   (ERR0R1)   1260,1280,1260 
1260  ERR0R1=ERR0R1 -ERROR 

IF   (ABS(ERRORl)-TOLER)  1290,1290,1270 
1270  SL0PE=DELTA1/ERR0R1 
1280  DELTA1=SL0PE*ERR0R 

ERR0R1=ERR0R 
1290  TEMP=TEMP+0ELTA1 

1300  IF   (PRESS-COEFFT(D)   1320,1320,1310 

1310   IF   ((.00001178*PRESS+.09411)*PRESS+382.1-TEMP)   1320,1320,1620 
1320  T=.55555556*TEMP+255. 38223 

TLOG=ALOG(T) 

TAU=1.0/T 

C0EFFT(222)=TAU*TAU*186210.06 

C0EFFT(217)=EXP(C0EFFT(222)+7.8791476-TL0G) 

C0EFFT(221)=162460.*TAU 

C0EFFT(220)=126970.*TAU 

C0EFFT(219)=EXP (149. 27765-23. *TLOG) 

C0EFFT(212)=1.89-C0EFFT(217) 

C0EFFT(211)=82.546-C0EFFT(221) 

C0EFFT(210)=.2I828*T-C0EFFT(220) 

COEFFT(209)=COEFFT(219)-.0003635*T 

C0EFFT(222)=(C0EFFT(222)+C0EFFT(222)+1.0)*C0EFFT(217) 

C0EFFT(217)=C0EFFT(212)-C0EFFT(222) 

C0EFFT(214)=C0EFFT(217)/C0EFFT(212) 

C0EFFT(216)=C0EFFT(211)*C0EFFT(214)+41.273-C0EFFT(221) 

C0EFFT(215)=C0EFFT(210)*C0EFFT(214)-.5*C0EFFT(220) 
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C0EFFT(214)=C0EFFT(209)*C0EFFT(214)+1.8461538*C0EFFT(219) 

COEFFT  (22 1)=.50*C0EFFT  (21 D-COEFFT  (216) 

C0EFFT(220)=.25*C0EFFT(210)-C0EFFT(215) 

COEFFT (219)  =  . 076923077*C0EFFT (209  )-COEFFT (214) 

C0EFFK2 18) =( .0037299965*T+14. 531813) *T+17806. 513+472. 2493*TL0G 

COEFFT (223 )=-.0074599931*T-14.531813*TL0G+2. 4524439+472. 24937 
**TAU 
1330  P=.068046190*PRESS 

PLOG=ALOG(PRESS) 

COEFFT (213)=C0EFFT(212)*TAU 

C0EFFT(226)=C0EFFT(213)*P 

COEFFT (225)=C0EFFT(226)*C0EFFT  (226) 

C0EFFT(224)=C0EFFT(225)*C0EFFT(225) 

COEFFT (224)=C0EFFT (224)*C0EFFT (224)*C0EFFT (226  ) 

C0EFFT(226)=C0EFFT(226)*C0EFFT(213) 

C0EFFT(213)=P*TAU 

STEAM(1)=(((C0EFFT(209)*C0EFFT(224)+C0EFFT(210))*C0EFFT(225)+ 
*C0EFFT(211))*C0EFFT(226)+C0EFFT(212)+4.55504/C0EFFT(213)) 
**. 0160185 

STEAM(2)=((((C0EFFT(214)*C0EFFT(224)+C0EFFT(215))*C0EFFT(225)+ 
*COEFFT (216) )*COEFFT (226 )+COEFFT (217) )*P+COEFFT (218))*. 04355685 

STEAM(3)=((((C0EFFT(219)*C0EFFT(224)+C0EFFT(220))*C0EFFT(225)+ 
*C0EFFT(221))*C0EFFT(226)+C0EFFT(222))*C0EFFT(213)+C0EFFT(223) 
*+4.55504*PL0G)*(-. 02419825) 

QUAL=1.0 

I=LINK+1 

GO  TO  (1640,1640,1340,1240,1340,1450),! 
1340  X=STEAM(II) 

IF  (H-X)  1350,1640,1360 
1350  Y=SATLIQ(II) 

QUAL=(H-Y)/(X-Y) 

X=1.0-QUAL 

STEAM(1)=QUAL*STEAM(1)+X*SATLIQ(1) 

STEAM(JJ)=QUAL*STEAM(JJ)+X*SATLIQ(JJ) 

STEAM(II)=H 

IF  (LINK-2)  1640,1640,1460 
1360  I=(II-1)*11+154 

J=l 

X=PLOG 
1370  DO  1380  K=209,219 

COEFFT(K)=COEFFT(I) 
1380       1=1+1 

X=X-C0EFFT(209) 

Y=Y-C0EFFT(210) 

C0EFFT(209)=(C0EFFT(211)*X+C0EFFT(212))*X+C0EFFT(213) 

COEFFT (210 )=(COEFFT (214) *X+COEFFT (215) )*X+COEFFT (216) 

C0EFFT(211)=( (C0EFFT(217)*X+C0EFFT(218) )*X+C0EFFT(219) )*Y 

Z=C0EFFT(211)+C0EFFT(210) 

SL0PE=C0EFFT(211)+Z 

Z=Z*Y+C0EFFT(209) 

GO  TO  (1390, 1430, 1530, 1540), J 
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1390  LINK=3 

IF  (Z-TEMP)  1300,1300,1400 
1400  TEMP=Z 

GO  TO  1300 
1410  J=2 

1=227 

X=H 

IF  (Z-H)  1420,1370,1370 
1420  1=187 

GO  TO  1370 
1430  DPDS=SLOPE 

PRESS=EXP(Z) 

PLOG=Z 

IF  (1-227)  1520,1520,1440 
1440  LIMK=4 

11=2 

JJ=3 

GO  TO  1120 
1450  DH=H-STEAM(2) 

STEAM(3)=ST£AM(3)+DH/(TEMP+459.69) 
1460  DS=S-STEAM(3) 

DELTA1=STEAM(3)-ERR0R1 

IF  (ERR0R1)  1470,1490,1470 
1470  IF  (ABS(DELTAl)-. 000005)  1500,1500,1480 
1480  DPDS=(PL0G-PL0G1)/DELTA1 
1490  PL0G1=PL0G 

ERR0R1=STEAM(3) 
1500  PRESS=PRESS*(1.0+DS*DPDS) 

IF  (LINK-4)  1510,1510,1550 
1510  IF  (ABS(DS)-. 000005)  1640,1120,1120 
1520  LINK=5 

X=H 

Y=S 

1=198 

J=3 

GO  TO  1370 
1530  TEMP=Z 

DTDS=SL0PE 

X=PL0G 

Y=H 

1=165 

J=4 

GO  TO  1370 
1540  DTDH=SL0PE 

GO  TO  1300 
1550  TEMP=TEMP+DS*DTDS 

IF   (ABS(DH)-.005)   1610,1560,1560 
1560  DELTA2=DH-ERR0R2 

IF   (ERR0R2)   1570,1590,1570 
1570  IF   (ABS(DELTA2)-.005)    1600,1600,1580 
1580  DTDH=(TEMP1-TEMP)/DELTA2 
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1590  TEMP1=TEMP 

ERR0R2=DH 
1600  TEMP=TEMP+DH*DTDH 

GO  TO  1320 
1610   IF   (ABS(DS)-.000005)   1640,1600,1600 
1620  WRITE   (6,1630)   NI,XI,YI 
1630  FORMAT   ( 13 ,2E14.8,13H  OUT  OF  RANGE) 

STOP 
1640  RETURN 

END 
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APPENDIX   B 
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BRAYTON  CYCLE:   PRESSURE  RATIO =2. 00 
WORKING  PRESSURE =135.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FLOW  (Lbm/s)=444.37 


STATE 
1 
2 

3 
4 

5 
6 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 
540.0 

675.7 
1668.8 
2060.0 
1804.2 

811.2 

784.0 

778.3 


DPP 
.078 


DHD 
.744 


SUMMARY: THERMAL   EFF =0.370 

HEATER  EFF=0.933  (PHI=0.90) 

SFC ( Lbm/HP-hr ) = 0 . 572 

VOL  T0TAL(-cu  ft)=5l6.5 

POWER (HP) =25000.0 

WORK  FRAC   BY  STM  CYCLE=0.045 


VOL(cu  ft) 

15.5      EFF=0.88 


.007      .260 

.0005 
.014 


10.0 

139.8 
29.3 


EFF(P)=0.91 
EFF(P)=0.88 

EFF=0.94 


HEATER    CYCLE:      PRESSURE  RATIO    =4.79 
TOTAL   PRESSURE   DR0P=.43 
AIR   FL0W(Lbm/s)=47.84 
FUEL  FLOW(Lbm/s)=    3-97 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1727.4 

4h 

5054.5 

5h 

1768.8 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP   DHD 

VOL(cu 

ft) 

RECUP 

.321    .650 

2.6 

EFF=0.96 

TURBINE 

1.1 

HEATER 

.012    .310 

67.4 

COMPRESSOR 

16.0 

EFF(P)=0.85 

RANKINE   CYCLE:      SATURATION   PRESSURE= 
STEAM  FLOW(Lbm/s)=      2.96 
FLOW  RATI0=0.0067 


56.7 


STATE 

TEMP(deg  R) 

PRESSURE (psi) 

Is 

540.0 

0.5 

2s 

555.0 

5^.7 

SF 

749.0 

56.7 

3s 

776.2 

56.7 

4s 

54o.o 

0.5  QUAL=0.844 

COMPONENT 

V0L(cu  ft) 

BOILER 

162.3 

TURBINE 

4.3  EFF(IS) 

=  0.88 

CONDENSER 

68.3 
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BRAYTON  CYCLE:   PRESSURE  RAT 10=  3.00 
WORKING  PRESSURE  =135.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FL0W(Lbm/s)=287.63 


STATE 
1 
2 

3 
4 

5 
6 

PP 

7 

COMPONENT 
REGEN 

TURBINE 

COOLER 
COMPRESSOR 

HEATER 

BOILER 


TEMP( deg  R) 

540.0 

769.8 
1538.8 
2060.0 
1643.6 

874.7 
841.2 

83I.8 


DPP 
.050 


DHD 
.837 


SUMMARY: THERMAL  EFF=0.432 

HEATER  EFF=0.933  (PHI=0.90) 

SFC(Lbm/HP-hr) =0.490 

VOL  TOTAL(cu  ft)=    318.1 

POWER (HP) =25000.0 

WORK  FRAC  BY  STM  CYCLE=0.044 


VOL(cu  ft) 

11.5     EFF=0.88 


.024      .352 

.0004 
.026 


5.3 

88.9 
20.9 


EFF(P)=0.91 
EFF(P)=0.88 
EFF=0.94 


HEATER    CYCLE:    PRESSURE   RATIO =  4.79 
TOTAL   PRESSURE   DR0P=.43 
AIR   FLOW(Lbm/s)=  40.99 
FUEL  FLOW ( Lbra/s ) =    3.40 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1586.5 

4h 

4932.9 

5h 

1638.8 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

VOL(cu 

ft) 

RECUP 

.310      .650 

2.2 

EFF=0.93 

TURBINE 

0.9 

EFF(P)=0.89 

HEATER 

.022      .358 

28.6 

COMPRESSOR 

12.7 

EFF(P)=0.85 

RANKINE   CYCLE:      SATURATION   PRESSURE=    128.0 
STEAM  FLOW ( Lbra/s ) =   2.48 
FLOW  RATI0=0.0086 


STATE 

TEMP(deg  R) 

PRESSURE (psi) 

Is 

540.0 

0.5 

2s 

555.0 

128.0 

SF 

806.2 

128.0 

3s 

839.7 

128.0 

4s 

540.0 

0.5     QUAL=0.819 

COMPONENT 

V0L(cu  ft) 

BOILER 

88.3 

TURBINE 

3.4      EFF(IS) 

=  0.88 

CONDENSER 

55-5 
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BRAYTON   CYCLE:      PRESSURE   RATIO =   4.00 
WORKING    PRESSURE   =    135-0 
TOTAL  PRESSURE  DR0P=.10 
MASS   FLOW    (Lbm/s)=2  39.91 


STATE 
1 
2 

3 

4 

5 
6 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP (d eg 

540.0 

843.9 

14-53.8 

2060.0 

1536.9 
927.0 
876.8 
860.6 

DPP 

.021 


R) 


S  UMMAR  Y : THERMAL   EFF =0.447 

HEATER  EFF=0.933  (PHI=0.90) 

SFC(Lbm/HP-hr)=0.473 

VOL  TOTAL(cu  ft)=  352.2 

POWER (  HP)  =25000.0 

WORK   FRAC   BY  STM   CYCLE=0.063 


DHD 
.930 


VOL(cu  ft) 


.032      .444 

.0003 

.04-6 


13.2 

3.6 
89.8 
18.4 


EFF=0.88 
EFF(P)=0.91 

EFF(P)=0.88 
EFF=0.94- 


HEATER    CYCLE:      PRESSURE  RATIO  =   4-. 79 
TOTAL   PRESSURE  DR0P=.43 
AIR  FLOW(Lbm/s)=39.59 
FUEL  FLOW(Lbm/s)=    3.29 


STATE 
lh 
2h 
3h 
4h 
5h 
6h 
7h 
COMPONENT 
RECUP 
TURBINE 
HEATER 
COMPRESSOR 


TEMP( d eg  R) 

535.0 

835.2 
14-94-.  5 
4855-0 
1553.8 

94-5-0 

760.O 
DPP        DHD 
.299       .650 

.033      .4-06 


V0L(cu  ft) 

2.1      EFF=0.92 
0.8     EFF(P)=0.89 

17.7 

12.0      EFF(P)=0.85 


RANKINE    CYCLE:      SATURATION   PRESSURE=    199-9 
STEAM  FLOW(Lbm/s)=      3. 19 
FLOW  RATI0=0.0133 


STATE 

TEMP(deg 

R) 

PRESSURE (psi) 

Is 

54-0.0 

0.5 

2s 

555.0 

199.9 

SF 

841.8 

199.9 

3s 

892.0 

199.9 

4-s 

540.0 

0.5   QUAL=0.810 

COMPONENT 

VOL(cu 

ft) 

BOILER 

119.0 

TURBINE 

5.1 

EFF 

CONDENSER 

70.5 
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BRAYTON  CYCLE:   PRESSURE  RATIO=  5.00 
WORKING  PRESSURES 35.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FLOW  (Lbm/s)=212.34 


STATE 
1 
2 


I 


5 
6 

PP 
7 
COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 

54o.o 

905.8 
1391.8 
2060.0 
1458.1 

972.1 

876.8 

846.9 
DPP        DHD 
.010      .950 

.022       .536 

.0003 
.068 


SUMMARY:      THERMAL   EFF=0.460 

HEATER   EFF=0.933    (PHI=0.90) 

SFC(Lbm/HP-hr)= 0.460 

VOL  T0TAL(cu  ft)=596.1 

POWER(HP) =25000.0 

WORK  FRAC   BY   STM  CYCLE=0.111 


V0L(cu  ft) 

15.7      EFF=0.88 


2.7 

107.8 

16.7 


EFF(P)=0.91 
EFF(P)=0.88 
EFF=0.94 


HEATER    CYCLE:       PRESSURE  RATIO =   4.79 
TOTAL  PRESSURE  DR0P=.43 
AIR   FL0W(Lbm/s)=    38.49 
FUEL  FLOW(Lbm/s)=    3-19 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1427.4 

4h 

4798.9 

5h 

1491.8 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

V0L(cu  ft) 

RECUP 

.289      .650 

2.1        EFF=0.90 

TURBINE 

0.8        EFF(P)=0.89 

HEATER 

.044      .454 

12.4 

COMPRESSOR 

11.5        EFF(P)=0.85 

RANKINE  CYCLE:   SATURATION  PRESSURE=  200.0 
STEAM  FLOW(Lbm/s)=  5-21 
FLOW  RATI0=0.0245 


STATE 

TEMP(deg  R) 

PRESSURE (psi) 

Is 

54o.o 

0.5 

2s 

555.0 

200.0 

SF 

841.8 

200.0 

3s 

937.1 

200.0 

4s 

54o.o 

0.5     QUAL=0.826 

COMPONENT 

VOL(cu  ft) 

BOILER 

298.1 

TURBINE 

10.9     EFF(IS)=0.88 

CONDENSER 

117.4 
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BRAYTON  CYCLE:   PRESSURE  RATIO  =6.00 
WORKING  PRESSURE =135.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FLOW (Lbm/s) =197.00 


STATE 
1 
2 

3 
4 

I 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 

540.0 

959-4 
13^3.8 
2060.0 
1396.2 
1011.8 

876.8 

835.3 


SUMMARY:  THERMAL   EFF= 0.464 

HEATER   EFF=0.933  (PHI=0.90) 

SFC(Lbm/HP-hr) =0.456 

VOL  TOTAL(cu  ft) =879-0 

POWER ( HP) =25000.0 

WORK  FRAC  BY  STM  CYCLE =0.151 


DPP 
.010 


DHD 
.950 


VOL(cu  ft) 


.016      .628 

.0003 
.07^ 


14.4 
2.2 

129.7 
15.8 


EFF=0.88 
EFF(P)=0.91 

EFF(P)=0.88 
EFF  =  0.9/+ 


HEATER  CYCLE:   PRESSURE  RATIO=  4.79 
TOTAL  PRESSURE  DR0P=.43 
AIR  FLOW (Lbm/s) =38. 17 
FUEL  FLOW (Lbm/s)-  3-17 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

4h 

1375.^ 

4755.8 

5h 

1443.8 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

VOL(cu  ft) 

RECUP 

.278      .650 

2.1        EFF=0.89 

TURBINE 

0.8        EFF(P)=0.89 

HEATER 

.055      .502 

9.5 

COMPRESSOR 

11.4        EFF(P)=0.85 

RANKINE   CYCLE:      SATURATION   PRESSURE=   200.0 
STEAM  FLOW ( Lbm/s )=   6. 70 
FLOW  RATIO =0.0 340 


STATE 

TEMP(deg  R) 

PRESS URE(psi) 

Is 

54o.o 

0.5 

2s 

555.0 

200,0 

SF 

841.8 

200.0 

3s 

976.8 

200.0 

4s 

540.0 

0.5     QUAL=0.838 

COMPONENT 

VOL(cu  ft) 

BOILER 

523.6 

TURBINE 

16.2     EFF(IS)=0.88 

CONDENSER 

153.3 
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BRAYTON  CYCLE:      PRESSURE  RATIO=7.00 
WORKING   PRESSURE =135.0 
TOTAL   PRESSURE  DRQP=.10 
MASS   Flow(Lbm/s) =187.48 


STATE 
1 
2 

3 

4 

5 
6 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 

540.0 
1007.0 
1304.9 
2060.0 
13^5.5 
10*47.6 

876.8 

825.2 


DPP 
.010 


DHD 
•  950 


SUMMARY: THERMAL  EFF=0.463 

HEATER   EFF=0.933  (PHI=0.90) 

SFC (Lbm/HP=hr) =0 .457 

VOL  T0TAL(cu  ft)=1195-7 

POWER (HP) =25000.0 

WORK  FRAC  BY  STM  CYCLE=0.189 


V0L(cu  ft) 

13.6  EFF=0.88 


.012   .720 

.0004 
.078 


1.9 

155.9 
15.3 


EFF(P)=0.91 
EFF(P)=0.88 
EFF=0.94 


HEATER  CYCLE:   PRESSURE  RATIO =  4.79 
TOTAL  PRESSURE  DR0P=.43 
AIR  FLOW(Lbm/s)=  38.22 
FUEL  FLOW(Lbm/s)=  3-17 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1333.3 

4h 

4721.2 

5h 

1404.9 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

V0L(cu  ft) 

RECUP 

.267      .650 

2.1        EFF=0.87 

TURBINE 

0.8        EFF(P)=0.89 

HEATER 

.066      .550 

7-7 

COMPRESSOR 

11.4        EFF(P)=0.85 

RANKINE  CYCLE:   SATURATION  PRESSURE=  200.0 
STEAM  FLOW(Lbm/s)=  7.92 
FLOW  RATI0=0.0422 


STATE 

TEMP(deg  R) 

PRESSURE (psi) 

Is 

540.0 

0.5 

2s 

555.0 

200.0 

SF 

841.8 

200.0 

3s 

1012.6 

200.0 

4s 

540.0 

0.5     QUAL=0.849 

COMPONENT 

V0L(cu  ft) 

BOILER 

782.4 

TURBINE 

21.0     EFF(IS)=0.88 

CONDENSER 

I83.6 
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BRAYTON  CYCLE:   PRESSURE  RAT  10=  8.00 
WORKING  PRESSURE=  135-0 
TOTAL  PRESSURE  DR0P=  .10 
MASS  FLOW(Lbm/s)  =  181.22 


STATE 
1 
2 

3 
4 

5 
6 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 
540.0 
1049.8 

1272.5 
2060.0 

1302.9 

1080.2 

876.8 

816. 3 
DPP   DHD 
.010  .950 

.010   .812 

.0004 
.080 


SUMMARY: THERMAL  EFF= 0.460 

HEATER  EFF=0.933  (PHI=0.90) 

SFC ( Lbm/HP-hr ) =  0 .459 

VOL  TOTAL(cu  ft)=1550.7 

POWER ( HP) =25000.0 

WORK  FRAC  BY  STM  CYCLS=0.225 


VOL(cu  ft) 

13.1      EFF=0,88 


1.6 
181.1 
15.1 


EFF(P)=0.91 
EFF(P)=0.88 

EFF=0.94 


HEATER   CYCLE:      PRESSURE  RAT  10=   4.79 
TOTAL   PRESSURE  DR0P=.43 
AIR  FLOW (Lbm/s) =38.45 
FUEL  FLOW(Lbm/s)=    3. 19 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1298.2 

4h 

4692.5 

5h 

1372.5 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

V0L(cu  ft) 

RE  CUP 

.256      .650 

2.1      EFF=0.86 

TURBINE 

0.8     EFF(P)=0.89 

HEATER 

.076      .598 

6.5 

COMPRESSOR 

1 

11.5     EFF(P)=0.85 

RANKINE   CYCLE:      SATURATION   PRESSURE=200.0 
STEAM  FLOW(Lbra/s)=   8.97 
FLOW  RATI0=0.0495 


STATE 

TEMP(deg  R) 

PRESSURE (psi) 

Is 

540.0 

0.5 

2s 

555.0 

200.0 

SF 

841.8 

200.0 

4s 

1045.2 

200.0 

540.0 

0.5     QUAL=0.859 

COMPONENT 

VOL(cu  ft) 

BOILER 

IO83.I 

TURBINE 

25.5     EFF(IS)=0.88 

CONDENSER 

210.2 
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BRAYTON  CYCLE:   PRESSURE  RATIO=  9-00 
WORKING  PRESSURE =135.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FLOW (Lbm/s)=l 76. 9 5 


STATE 
1 
2 


I 


I 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 

540.0 
IO89.O 
1245.0 
2060.0 
1266.2 
1110.2 

876.8 

808.2 


SUMMARY: THERMAL  EFF=0.456 

HEATER  SFF=0.933  (PHI=0.90) 

SFC(Lbm/HP-hr)= 0.464 

VOL  TOTAL(cu  ft)=1950.3 

POWER (HP) =25000.0 

WORK  FRAC  BY  STM  CYCLE =0.260 


DPP 
.010 


DHD 
.950 


.010  .904 

.0004 
.080 


VOL(cu  ft) 
12.7 
1.4 

197.5 
15.0 


EFF=0.88 
EFF(P)=0.91 

EFF(P)=0.88 
EFF=0.94 


HEATER  CYCLE:   PRESSURE  RATIO=  4.79 
TOTAL  PRESSURE  DROP  =  .4 3 
AIR  FLOW(Lbm/s)=  38. 80 
FUEL  FLOW(Lbm/s)=  3-22 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1268.3 

4h 

4668.2 

5h 

1345.0 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

V0L(cu  ft) 

RECUP 

.245      .650 

2.2      EFF=0.85 

TURBINE 

0.8      EFF(P)=0.89 

HEATER 

.087      .646 

5.7 

COMPRESSOR 

11.7      EFF(P)=0.85 

RANKINE  CYCLE:   SATURATION  PRESSURE=  200.00 
STEAM  FLOW(Lbm/s)=  9.90 
FLOW  RATIO=0.0560 

STATE  TEMP(deg  R)  PRESSURE ( psi) 

Is  540.0  0.5 

2s  555.0  200.0 

SF  841.8  200.0 

3s  1075.2  200.0 

4s  540.0  0.5  QUAL=0.867 

COMPONENT  V0L(cu  ft) 

BOILER  1439.0 

TURBINE  29.8  EFF(IS)  =0.88 

CONDENSER  234.5 
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BRAYTON  CYCLE:   PRESSURE  RAT 10 =10. 00 
WORKING  PRESSURE =135.0 
TOTAL  PRESSURE  DR0P=.10 
MASS  FLOW  (Lbm/s)=174.00 


STATE 
1 
2 


I 


5 
6 

PP 

7 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


TEMP(deg  R) 
540.0 
1125.1 
1221.1 
2060.0 
1234.2 

1138.1 
876.7 

800.9 
DPP        DHD 
.010     .950 

.008      .996 

.0004 
.082 


SUMMARY: THERMAL     EFF=0.451 

HEATER   EFF=0.933   (PHI=0.90) 

SFC ( Lbm/HP-hr ) =0 . 469 

VOL  TOTAL(cu  ft) =2403.0 

POWER (HP) =25000.0 

WORK  FRAC  BY  STM  CYCLE=0.296 


VOL(cu  ft) 

12.4     EFF=0.88 


1.3 

232.8 

14.9 


EFF(P)=0.91 
EFF(P)=0.88 

EFF=0.9^ 


HEATER   CYCLE:      PRESSURE  RATI0=   4.79 
TOTAL   PRESSURE  DR0P=.43 
AIR   FL0W(Lbm/s)=    39.22 
FUEL  FLOW(Lbm/s)=    3.25 


STATE 

TEMP(deg  R) 

lh 

535.0 

2h 

835.2 

3h 

1242.5 

4h 

4647.3 

5h 

1321.1 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP        DHD 

V0L(cu  ft) 

RECUP 

.235      .650 

2.2     EFF=0.84 

TURBINE 

0.8     EFF(P)=0.89 

HEATER 

.098      .694 

5.0 

COMPRESSOR 

11.9      EFF(P)=0.85 

RANKINE   CYCLE:      SATURATION   PRESSURE=    199-9 
STEAM  FLOW(Lbra/s)=    10. 76 
FLOW  RATI0=0.0618 


STATE 

TEMP(deg 

R) 

PRESSURE (psi) 

Is 

540.0 

0.5 

2s 

555.0 

199.9 

SF 

841.7 

199.9 

3s 

1103.1 

199.9 

4s 

540.0 

0.5     QUAL=0.875 

COMPONENT 

V0L(cu 

ft) 

BOILER 

I830.6 

TURBINE 

34.0 

EFF( 

=  0.88 

CONDENSER 

257.1 
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